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SECTION  1  -  SUMMARY 


This  grant  covers  the  Air  Force  Research  Initiation  Grant  (AFERIG)  Program  for  the  Academic 
Year  1990-1991.  It  was  announced  in  the  Fall  of  1990  for  up  to  25  AFERIGs  carrying  a  stipend 
of  $20,000.  The  Engineering  Foundation’s  five  Founder  Societies  received  a  total  of  262 
proposals.  Fifty-six  of  these  were  selected  for  referral  to  AFOSR.  AFOSR  selected  25  with  10 
alternates.  Eight  of  the  25  also  received  various  NSF  grants  and  were  replaced  by  alternates. 
The  final  list  of  awardees  is  found  in  Section  2.  All  awardees  have  submitted  reports  for  this  final 
technical  report.  Abstracts  of  these  studies  are  included  in  Section  3. 
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Ultrasonic-Based  Measurement  of  Residual  Stress  Induced  by  Large 

Deformation  Processes 


Professor  Jahan  Rasty 

Department  of  Mechanical  Engineering 
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ABSTRACT 


Evaluation  of  the  residual  stress  in  a  material  is  one  of  the  fundamental  tasks 
of  a  mechanical  engineer.  Conventional  methods  of  stress  determination  include 
hole  drilling,  sectioning,  and  X-ray  scanning.  A  major  drawback  of  these  methods 
is  that  they  are  either  destructive  or  restricted  to  surface  measurement  only. 
Efficient  utilization  of  high-tech  alloys  in  complex  metal  forming  processes  requires 
development  of  simple  techniques  for  determining  the  state  of  residual  stress.  The 
application  of  ultrasonic  waves  to  material  testing  and  characterization  is  a 
promising  approach  that  does  not  have  the  aforementioned  drawbacks. 

Specifically,  ultrasonic  technique  has  the  advantage  of  being  nonhazardous, 
nondestructive,  versatile,  and  applicable  to  through  thickness  measurements. 

The  objective  of  this  research  was  to  further  develop  the  acoustoelasticity 

theory  for  the  application  to  the  residual  stress  analysis.  In  the  first  phase  of  this 

work,  a  comprehensive  literature  search  was  conducted  to  determine  the  state  of 

technology  in  this  field.  Although  a  considerable  amount  of  research  has  been 

done  in  this  area,  most  of  it  was  scattered.  This  research  presented  a  unified 

treatment  of  the  acoustoelasticity  theory  for  application  to  residual  stress  analysis. 

Furthermore,  a  modification  to  the  small  plastic  deformation  theory  was  proposed 

•* 

for  application  to  large  plastic  deformation  as  encountered  in  manufacturing 
processes.  An  ultrasonic  experimental  setup  was  developed  for  residual  stress 
analysis  utilizing  the  pulse-echo-overlap  technique.  Further  experiments  were 
outlined  for  determining  the  validity  of  the  proposed  theory. 
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RI-B-90-02 


The  Effects  of  Grain  Size  on  Cyclic  Deformation  and  Fatigue 

Damage  Accumulation 


Professor  David  J.  Morrison 

Department  of  Mechanical  and  Aeronautical  Engineering 
Clarkson  University 
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Executive  Summary 


This  report  describes  work  accomplished  during  a  one  year  Air  Force  Engineering 
Foundation  Grant,  RI-B-90-02,  to  study  the  effects  of  grain  size  on  cyclic  deformation  and 
fatigue  damage  accumulation.  Significant  results  have  been  obtained  in  the  project. 
Accomplishments  include  the  development  of  a  computer  program  that  automatically 
calculates  key  parameters  describing  cyclic  stress-strain  behavior  and  fundamental 
mechanisms  of  cyclic  plasticity,  a  correlation  between  hysteresis  loop  shape  behavior  and 
the  evolution  of  surface  fatigue  damage,  image  analysis  experiments  of  surface  damage 
evolution,  crystallographic  texture  experiments,  and  a  correlation  between  strain 
localization  phenomena  and  fatigue  crack  initiation. 

The  results  of  the  studies  indicate  that  cyclic  strain  localization  phenomena  can  be 

accurately  measured  by  evaluating  subtle  changes  in  the  shapes  of  the  cyclic  stress-strain 

hysteresis  loops.  Cyclic  strain  localization  was  found  to  be  much  more  pronounced  at  low 

plastic  strain  amplitudes  where  persistent  slip  bands  (psbs)  form.  Furthermore,  the  strain 

localization  is  more  pronounced  as  the  grain  size  of  the  material  increases.  These 

* 

observations  of  strain  localization  phenomena  have  significant  engineering  significance. 
Specifically,  at  low  plastic  strain  amplitudes,  an  increase  in  grain  size  has  a  detrimental 
effect  on  fatigue  life.  At  high  plastic  strain  amplitudes,  where  plastic  strain  is  not  localized, 
grain  size  has  very  little  effect  on  fatigue  life.  These  research  findings  have  resulted  in  one 
journal  publication  and  another  manuscript  in  preparation  for  publication.  A  reprint  of  the 
journal  article  has  been  included  in  Appendix  111. 
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RI-B-90-03 


Phase  Transformations  and  Microstructure  Development  in  Two-Phase 
TiAl  +  Ti3Al  Zntermetallic  Alloys 


Professor  Vijay  Vasudevan 

Department  of  Materials  Science  and  Engineering 
University  of  Cincinnati 
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COMPOSITION  DEPENDENCE  OF  THE  MASSIVE  TRANSFORMATION 
FROM  a  TO  y  IN  QUENCHED  TiAl  ALLOYS 


Ping  Wang  and  Vijay  K.Vasudevan 

Dept,  of  Materials  Science  and  Engineering,  University  of  Cincinnati,  ML  12,  Cincinnati,  OH  45221 

(Received  May  5,  1992) 


Introduction 


Alloys  based  on  the  intermetallic  compound  TiAl  have  been  of  much  interest  in  recent  years  as  lightweight, 
structural  materials  for  elevated  temperature  aerospace  applications  (1).  The  most  promising  alloys  arc  based  on  the  Ti-48 
A1  composition  (compositions  in  the  text  refer  to  at.%)  with  ternary  or  quartemary  additions  and  with  the  two-phase 
TiyAI+TiAl  (a2+y)  lamellar  microstructure  (2-5).  Several  studies  concerned  with  phase  equilibria  and  transformations  in 
Ti->40  A1  alloys  have  been  reported  (6-10).  In  low  A1  alloys,  the  lamellar  structure  has  been  suggested  to  be  produced  by 

the  reaction  a-»oc2-»a2+y,  involving  initial  formation  of  stacking  faults  (11).  More  recently,  it  has  been  suggested  that 
this  mechanism  is  operative  till  -42  Al,  whereas  at  higher  concentrations  the  y  lamellae  are  formed  directly  from  the  a  and 
subsequently  the  a  orders  to  aj  to  produce  the  same  morphology  of  lamellae  as  observed  at  lower  concentrations  (4,12). 
Very  recently,  in  two  studies  it  has  been  shown  that  cooling  rate  has  a  major  effect  on  the  transformation  of  the  a  phase. 
In  the  First,  alloys  containing  47-48  Al  were  studied  (13),  and  it  was  observed  that  the  y  forms  from  the  a  in  three 
morphologies  depending  on  cooling  rate.  At  low  cooling  rates,  the  lamellar  morphology  prevails,  at  intermediate  rates  a 
Widmanstatten  morphology  appears  and  at  very  high  cooling  rates  a  massive  transformation  dominates.  The  occurrence  of 
the  latter  in  a  quenched  Ti-48  Al  alloy  has  also  been  confirmed  very  recently  (12).  In  the  second,  where  a  Ti-48  Al-2Mn- 
2Nb  alloy  was  studied  (14),  very  similar  effects  of  cooling  rate  on  the  microstructure  were  reported,  and  the  occurrence  of 
the  massive  reaction  at  high  cooling  rates  was  confirmed  (15).  The  purpose  of  this  paper  is  to  report  on  the  composition 
dependence  of  transformations  in  Ti-(43-50)  Al  alloys  water  quenched  from  the  fully  a  region.  In  particular,  evidence  for 
a  massive  transformation  from  a  to  y  in  alloys  with  46.54  to  50  Al  is  presented. 

Experimental 

The  compositions  of  the  alloys  studied  are  shown  in  Table  1.  The  Ti-43,  46  and  50  Al  alloys  were  prepared  by 
arc-melting  commercial  purity  Ti  and  99.999%  pure  Al.  The  48  Al  alloy  was  obtained  as  ingots  in  the  cast+HIP  (1200°C, 
2h,  105  MPa)  condition  from  commercial  sources.  The  43,  46  and  50  Al  alloys  were  homogenized  at  1 100°C  for  50  h 
followed  by  slow  cooling.  Cut  samples  from  the  alloys  were  coated  with  a  high  temperature  glass  ceramic  (Coming  Glass 
7913)  for  protection  from  the  atmosphere  and  heated  to  the  completely  a  region,  namely,  1325, 1360, 1400°C  and  1435- 
1450°C  for  the  43,  46,  48  and  50  Al  alloys,  respectively.  After  holding  for  various  times  from  15  minutes  to  1  h,  the 
samples  were  cooled  at  different  rates,  namely,  furnace  cool  (FC),  air  cool  (AC),  oil  quench  (OQ),  water  quench  (WQ) 
and  iced  brine  quench  (IBQ).  Analysis  of  the  oxygen  and  nitrogen  contents  after  these  various  heat  treatments  showed  that 
any  pick-up  of  oxygen  and  nitrogen  was  restricted  to  the  outer  surface  of  the  samples,  whereas  the  interiors  were 
unaffected  (Table  1  and  ref.  (13)).  Analysis  of  the  Al  contents  after  heat  treatment  showed  that  a  small  loss,  typically 
<0.25  at.%,  had  occurred  in  the  alloys. 


TABLE  1 

Chemical  Compositions  of  Alloys  Studied 


Aimed  Composition  (at.%) 

Actual  Composition 

Ti  (at.%) 

Al  (at.%) 

O  (wt.%) 

N  (wt.%) 

Ti-43A1 

0.062 

■oisb 

Ti-46AI 

1 

0.051 

Ti-46A1  HT* 

I 

1BK  ilHiKfe-il. 

0.052 

0.005 

Ti-48  Al 

■ 

47.86 

0.049 

0.006 

Ti-48  Al  HT* 

1 

47.98 

0.036 

0.004 

Ti-50A1 

Bal 

49.40 

0.088 

0.014 

*Post  heat  treatment  analysis  after  removal  of  -2mm  of  outer  surface  of  sample 
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The  microstructures  of  the  alloys  were  observed  optically  am!  also  in  a  JKOL  JXA-840  scanning  electron 
microscope  (SEM)  utilizing  back  scattered  electrons  (USE)  for  obtaining  atomic  number  contrast.  Fully  quantitative 
electron  microprobe  analysis  of  the  compositions  of  phases  was  conducted  in  the  same  SEM.  Thin  foils  for  transmission 
electron  microscopy  (TEM)  were  prepared  from  the  samples  using  conditions  described  previously  (11).  'lire  foils  were 
observed  in  a  Philips  CM20  TEM  operated  at  200  kV,  utilizing  bright  field  (HE),  dark  field  (HE)  and  selected  area 
diffraction  (SAD)  imaging  modes. 


Results 

In  this  paper,  only  the  results  of  microstructures  of  WQ  samples  are  reported.  Thus,  optical  micrographs  of 
samples  of  the  various  alloys  WQ  from  the  completely  a  region  are  shown  in  the  micrographs  in  Figs.  la-d.  The  optical 
micrograph  of  the  WQ  Ti-43.19  A1  ahoy  appears  featureless.  Fig.  la,  and  TEM  examination  showed  that  the  micro¬ 
structure  was  composed  fully  of  02  domains.  A  prominent  feature  seen  in  the  micrograph  of  the  Ti-46.54  Ai  alloy.  Fig. 
lb,  is  the  presence  of  dark,  deeply  etched,  patchy  regions  in  a  light,  featureless  matrix.  These  regions  appear  to  have 
nucleated  largely  at  the  prior  a  grain  boundaries  and  then  spread  rapidly  info  the  grains  in  a  spheroidal  and  massive 
fashion.  On  occasion,  these  dark  regions  were  also  found  within  grains,  at  some  distance  from  the  grain  boundaries. 
Closer  inspection  of  the  latter  sometimes  reveals  the  existence  of  a  Tine  striated  structure  emanating  from  the  grain 
boundaries  and  terminating  at  locations  where  the  dark  regions  commence.  The  interface  between  the  dark  and  light 
regions  varies  from  smooth  to  jagged  from  one  region  to  another.  The  exact  nature  of  the  interface,  i.e.,  incoherent  or 
semi-coherent,  has  not  been  determined  yet.  Examination  of  the  47.86  and  49.40  Al  alloys.  Figs,  lc-d,  revealed  that 
practically  the  entire  microstructure  of  the  samples  was  dominated  by  the  dark,  mottled  regions  seen  in  the  46.54  AI  alloy. 
A  feathery  morphology  was  also  present  often.  It  should  be  noted,  however,  that  at  the  higher  Al  contents,  particularly 
49.40  Al,  it  was  difficult  to  suppress  the  formation  of  the  "equilibrium"  lamellar  morphology,  which  was  often  present  at 
the  grain  boundaries  and  impinged  on  the  dark  regions  within  the  grains. 

TEM  observations  establish  that  the  light,  featureless  regions  seen  in  the  optical  micrograph  of  the  WQ  46.54  AJ 
alloy  consist  of  02  domains.  Fig.  2a,  whereas  the  dark  regions  are  y,  as  shown  in  Fig.  2b  by  the  BF  micrograph  and  inset 
H.IO]  y  SAD  pattern.  TEM  observations  of  the  47.86  and  49.40  Al  alloys,  revealed  that  practically  the  entire 
microstructure  was  dominated  by  y,  although  a  few,  sparsely  distributed  regions  of  lamellar  oq+y  were  observed.  An 
example  is  shown  in  Fig.  3  for  the  47.86  Al  alloy.  The  massive  y  phase  in  all  of  these  WQ  alloys,  especially  the  46.54 
Al,  is  typically  characterized  by  a  complex  substructure  consisting  of  antiphase  boundaries  (APB’s),  stacking  faults, 
microtwins  and  dislocations  (Figs.  2b  and  3).  The  presence  of  APB’s  might  suggest  the  occurrence  of  an  intermediate 
disordered  f.c.c.  or  f.c.t.  structure  prior  to  the  formation  of  y.  Electron  diffraction  analysis  of  many  y  grains  showed  the 
occurrence  of  all  six  variants.  Aside  from  f.c.c.-twin  related  variants,  numerous  contacting  variants  in  perpendicular  twin 
relation,  i.e.,  mutually  perpendicular  c-axes,  are  observed.  These  observations  suggest  that  a  direct  transformation  from  a 
to  y  had  occurred,  perhaps  in  a  massive  fashion. 

BSE  images  of  the  WQ  46.54  Al  sample  (Fig.  4a)  are  similar  to  those  reported  previously  in  the  47.20  Al  alloy 
(13)  and  show  regions  of  both  darker  contrast  corresponding  to  y  and  slightly  lighter  contrast  corresponding  to  aj- 
Although  the  observation  of  contrast  differences  between  the  two  phases  indicates  a  difference  in  chemical  composition,  a 
microprobe  scan  analysis  of  the  02  and  y  regions  extending  through  the  interface,  over  a  distance  of  -600  pm  at  -10  pm 
intervals.  Fig.  4b,  reveals  that  the  compositions  of  the  two  phases  are  practically  identical.  Furthermore,  no  consistent 
trend  in  the  compositions  of  the  two  phases  is  seen  that  might  explain  the  differences  in  contrast  seen  in  the  BSE  images. 
The  possible  reasons  for  the  origin  of  the  contrast  differences  have  been  discussed  elsewhere  (13).  Observations  of  the 
47.86  Al,  and  49.40  Al  alloys  using  BSE  imaging  did  not  reveal  any  contrast,  i.e.,  were  uniformly  dark  and  no  interfaces 
(between  any  a2  and  y)  could  be  discerned.  Microprobe  scan  analysis  of  these.  Figs.  5a-b,  respectively,  shows  that  the 
compositions  are  practically  invariant  and  near  that  of  the  bulk  alloy  compositions.  These  observations,  together  with  the 
TEM  results,  provide  strong  evidence  for  a  massive  transformation  from  a  to  y  at  high  cooling  rates. 

Discussion 


As  mentioned  earlier,  in  a  previous  study  (13)  the  effect  of  cooling  rate  (and  therefore,  undercooling)  on  the 
transformation  from  a  to  7  in  a  cast  and  homogenized  Ti-47.20  Al  alloy  and  a  cast  and  extruded  Ti-47.86  Al  alloy  were 
reported.  In  that  study  it  was  found  that  at  low  cooling  rates  the  7  phase  forms  in  a  lamellar  morphology,  at  intermediate 
rates  a  Widmanstatten  morphology  appears  and  at  very  high  rates  a  massive  morphology  is  prevalent.  In  the  present 
study,  the  effect  of  Al  content  in  binary  alloys  has  been  examined  and  the  results  show  that  the  massive  reaction  from  a  to 
7  occurs  at  high  cooling  rates  from  46.54  to  50  Al.  The  massive  transformation  has  been  extensively  studied  in  the  past 
and  examples  of  this  transformation  can  be  found  in  many  systems  (16-19).  Essentially,  a  massive  transformation  takes 
place  by  diffusional  nucleation  and  growth  and  involves  a  change  in  crystal  structure  without  a  change  in  composition, 
and  in  principle,  can  be  observed  in  any  system  where  the  pertinent  phase  fields  overlap  in  composition  (16,17). 
Moreover,  it  is  possible  for  a  massive  reaction  to  be  considered  as  occuring  in  the  metastable  two-phase  field,  such  as  that 
below  a  invariant  reaction  temperature  (eg.,  eutectoid  temperature)  (17).  Both  the  TEM  observations  showing  the 
structural  change  and  the  microprobe  analysis  establishing  composition  invariance  attest  to  the  occurrence  of  the  massive 
reaction  from  a  to  7  at  high  cooling  rates.  The  nucleation  of  the  massive  7  is  seen  to  occur  mostly  at  the  grain  boundaries 


12 


i 


m 


[;K i  I  Optical  tnic  rostmctmes  nl  1  iAI  alloys  WO  liom  ilic  lulls  <>  rep  ion  (a )  1  i  -1 '  1 1 1  A I  1 1' 
Al.  ami  <d>  I'i  4‘)  40  Al 


I  i  40  SI  Al.  (c)  1 1  4  S 


of  the  prior  a.  as  is  seen  usually  m  mosi  massive  reactions  ( I  h  I  *0  I  urtheniiore.  the  present  ohsei  v  at  ions  <>l  I  lie  masse 
reaction  tor  Al  contents  trout  SO  to  40  S  I  Al  imlu  ates  that  sipmticain  pen--tiation  into  the  equilibrium  toym  inetasiahl 
equilibrium  my  restion  occurs  However.  aitliouph  all  results  have  not  been  presented  lieiem.  alloys  eontaimne  1 '  I-  4 
WO  front  the  u  ret!  ion  1 1  ausf  oi  met  I  romp  letch  to  u  >  and  show  ci)  no  ev  idence  ot  the  masse,  e  1 1  ansi  ot  mat  ion 

In  ,  inlet  to  esp'.un  the  ob.etved  i  omposinon  dependeiiee  "l  the  massive  tea.  nmi.  ai »  mate  Iiio'a!  , 

eqmrihtmm  ami  melastahlc  eqmlihiuim  phase  houndatie-s  and  teinpeiatuie  ot  mu  leatmn  d  the  in  tssive  ale  i  ; 

Ih  us .  the  tele  vant  port  ion  ot  cm  tent  I  \  accepted  I  i  A I  phase  i  hue  t  am  1 1  on  i  the  w  <  <t  k  ot  M,  t  u  I  loin'll  el  d  i 1 1  t,e . 

shehth  h\  I  ones  and  h  ant  maun  1  I  I .  is  show  n  in  lie  h  I  lie  esl  i  mated  pos  it  ions  ot  the  met  ust  able  d, .  ■  i  •  .  1 1 

phase  b  vmd.uies  and  the  n,d  and  "/<’  •  I  .  line  s  a :  e  a  Iso  shove  u  It  i.  an  bt*  seen  I  tia!  t  ot  a.  ■ .  o\  -\1  v  ontents  ap|c  ' '  ; 

■Jo  at  o  the  i  on  i  pos  1 1.  oi;s  ot  die  >i  and  phases  oveilap  above  the  elite  >  toid  !•  mi  petal  <  I  I  ,’s  <  1.  whereas 

i  i  an  ent  i  at  ions  i  he  oveilap  ov  cm  s  hrlow  the  ei  ik\  ton!  uunpei  at  me  in  l  he  me  last  a  hie  i  •  -  phase  field  1,  e, 
believed  that  the  I.,  line  lepiesenls  ll;e  limit  tot  the  massive  reaction  and  the  dtsiume  between  lb’-  I,,  hue  an  !  d 
tempetature  of  start  ot  die  reaction  represents  the  linden  nnhiip  and  tbendoie  die  o.nv  me  lon  e  loi  the  tn.issi ■. 
ti  at  is  toiuu!  i  on  <  In.  I  i  1 1 1  lie  tt  (  .’( 1 1  has.  Iiowevel.  at  mu  d  that  the  m.c  eve  tr.iu  latma'ia::  leqm’i's  linden  online'  '■  no 

,  I :  ;  1  *  ...  O,  ,•  |  ,  \i  alio-,  it-  ma  e  IV  e  lea:  I  ion  I  111  o.  m:t  h\  under:  nohn"  VV  U 


iov  \1  contents  a 

:■■  '  I  1  ,’s  Cl.  w  tie 
i  •  •■  phase  field  I 
e  be!  w een  dm  I  , 
it  iv  me  lotee  toi 
n  tequoes  linden 
,  m  h\  under:  no! 


■n  aaliii:-  "■  no 
mini "  •'  vv  d 


riu  oi  mu'  t 1 


COMPOSITION  (AT  PCT  Al) 


Vo 1 .  27,  No.  1 


MASS  1  VI-;  T RANS FORMAT  1  ON 


FIG.  4.  (a)  BSE  micrograph  and  (b)  microprobe  scan  analysis  of  Ti-46.54  Al  alloy  WQ  from  the  fully  a  region. 
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FIG.  5.  Microprobe  scan  analysis  of  (a)  Ti-47.86A1  and  (b)  Ti-49.40  alloys  WQ  from  the  folly  a  region. 


FIG  6.  Partial  phase  diagram  of  the  Ti-AI  system  (from  ref  (12)). 
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ssuming  for  the  sake  of  simplicity  that  the  a  to  ai  transformation  temperature  does  not  change  markedly  with  cooling 
ate,  the  temperature  range,  AT,  available  for  nucleation  and  growth  of  the  massive  reaction  from  a  to  yean  be  estimated 
rom  the  difference  between  T0  ajy  and  T0  a/a2.  Thus,  the  values  of  AT,  which  increase  rapidly  with  A1  content,  are 
stimated  to  be  approximately  125,  200,  265  and  335°C  for  the  45,  46,  47  and  48  Al  alloys,  respectively.  For 
ompositions  below  -44  Al,  the  aja 2  boundary  is  encountered  prior  to  the  ajy  T0  line  and  so  rapid  quenching  results  only 
(i  the  formation  of  02.  The  absence  of  the  massive  reaction  between  46  and  -44  Al  is  most  likely  related  to  either  the 
lifficulty  in  nucleating  y  at  these  compositions  or  that  the  available  AT  does  not  provide  an  adequate  driving  force  for  its 
tucleation  (before  the  a  to  a2  transformation  intervenes). 

Another  interesting  aspect  seen  is  that  the  microstructures  of  the  46.54  and  47.20  Al  (13)  alloys  consists  of  012  and 
nassive  y,  i.e.,  the  massive  transformation  does  not  appear  to  proceed  to  completion.  In  comparison  the  47.86  Al  alloy 
nitially  in  the  cast  +  HIP  condition  (present  study)  as  well  as  in  the  extruded  condition  (13)  appears  to  transform  almost 
:ompletely  to  massive  y.  It  was  thought  in  the  previous  study  that  this  difference  in  behavior  maybe  due  to  a  grain  size 
:ffect,  i.e.,  the  47.20  Al  alloy  had  a  prior  grain  size  >500  gm  compared  with  40  to  100  jim  in  the  extruded  47.86  Al  alloy 
ind  as  a  consequence  there  were  fewer  sites  for  nucleation  of  the  massive  y.  Once  nucleated  these  grow  until  the  0/02 
>oundary  is  encountered  below  1125°C  and  the  a-Ki2  transformation  intervenes  and  arrests  further  progression  of  the 
nassive  reaction.  For  the  reasons  discussed  in  the  preceding  paragraph,  it  is  clear  that  the  difference  in  alloy  compositions 
:ou!d  also  be  responsible  for  this  effect,  i.e.,  the  available  AT  (and  therefore  time)  for  nucleation  and  growth  of  the 
nassive  reaction  before  the  intervention  of  the  a  to  a2  transformation  decreases  markedly  with  a  decrease  in  the  Al 
tontent.  The  relative  contributions  of  these  two  effects  needs  investigation.  Finally,  it  is  clear  that  before  the  temperature 
lependence  of  the  massive  reaction  can  be  ascertained  positively,  experimental  determination  of  the  temperature  of  start  of 
he  reaction  is  required.  These  aspects,  together  with  the  nature  of  the  parent-massive  interface,  orientation  relationships, 
tnd  the  nucleation  and  growth  mechanisms  of  the  different  morphologies  of  y  are  the  subject  of  ongoing  research. 

Conclusions 

In  summary,  the  microstructures  of  Ti-(43-50)  Al  alloys  water  quenched  from  the  fully  a  phase  region  have  been 
studied.  In  particular,  evidence  for  the  occurrence  of  a  massive  reaction  from  a  to  y  in  alloys  with  Al  contents  in  the  range 
of  46.54  to  50  at.%  has  been  presented  and  rationalized  in  terms  of  phase  diagram  considerations. 
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Observation  of  a  Massive 
Transformation  from  a  to  y  in 
Quenched  Ti-48  At.  Pet  Al  Alloys 

PING  WANG,  G.B.  VISWANATHAN, 
and  VUAY  K.  VASUDEVAN 

Owing  to  their  unique  combination  of  attractive  prop¬ 
erties,  such  as  lightness,  superior  strength,  creep  and  ox¬ 
idation  resistance,  alloys  based  on  the  intermetallic 
compound  TiAl  have  been  of  much  interest  in  recent  years 
as  structural  materials  for  elevated  temperature  aero¬ 
space  applications.1'*  The  alloys  developed  to  date  that 
show  the  most  promise  with  regard  to  room-temperature 
ductility  and  elevated  temperature  strength  are  based  on 
the  Ti-48A1  composition  (compositions  in  the  text  refer 
to  atomic  percent),  with  additions  of  elements  such  as 
V,  Cr,  Nb,  Mn,  etc.,  singly  or  in  combination.*2-6*  All 
of  these  alloys  are  two-phase,  with  a  microstructurc 
composed  of  either  fully  lamellar  grains  of  TijAl  +  TiAl 
(or2  +  y)  or  duplex  primary  y  plus  lamellar  grains,  de¬ 
pending  on  processing  and  heat  treatment.*4*  The  prop¬ 
erties  of  these  alloys  are  quite  sensitive  to  microstructure 
and,  therefore,  to  heat  treatment;  duplex  structures  show 
higher  tensile  ductilities,  .for  example,  up  to  2.2  pet  in 
binary  Ti-48A1P-4*  and  up  to  4  pet  in  ternary  or  quarter- 
nary  alloys,**'6'  whereas  fully  lamellar  structures  have 
poor  ductility.*4-7' 

The  lamellar  (a2  +  y)  structure  was  first  studied  by 
Blackburn,'8*  who  reported  the  formation  of  alternate  la¬ 
mellae  of  y  and  a2  with  an  orientation  relationship  in¬ 
volving  {111}  y/{0001}  a2  and  (Oil)  y/(Tl20)  a2  and 
suggested  that  the  y  may  form  by  a  stacking  fault  mech¬ 
anism.  More  recently,  several  studies  of  phase  equilibria 
and  transformations  in  two-phase  TijAl  +  TiAl  alloys 
have  appeared.,,-,7'  Of  particular  significance  are  the  re¬ 
sults  of  these  investigations  showirfg  decidedly  that  the 
(/3  +  y)  two-phase  field  shown  in  previous  Ti-Al  phase 
diagrams"81  does  not  exist  and,  instead,  the  (a  +  y)  phase 
field  extends  up  to  the  peritectic  temperature;  i.e.,  the  a 
region  exists  up  to  1450  °C  and  50  at.  pet  Al.  Recently, 
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data  of  Antony151  to  account  for  all  of  the  observed  inter- 
planar  spacings  and  hence  the  recharacterization  of  the 
silicide  in  that  alloy. 

The  calculated  interplanar  spacings  for  the  silicides  5 , . 
and  S'  are  compared  with  the  values  determined  ex¬ 
perimentally  (X-ray  data  of  Antony151)  in  Table  I.  Now 
when  the  measured  interplanar  spacings  for  5'  silicide151 
are  compared  with  the  calculated  values,  it  may  be  noted 
that  the  experimentally  observed  interplanar  spacing  of 
0.350  nm  is  considerably  different  from  the  correspond¬ 
ing  calculated  value  of  0.335  nm.  Further,  the  inter¬ 
planar  spacings  0.303  and  0.201  nm  are  not  accounted 
for  by  either  S'  silicide  or  a-titanium.  However,  it  may 
be  noted  that  the  interplanar  spacings  0.350,  0.303,  and 
0.201  nm,  which  do  not  agree  for  S',  are  the  typical 
interplanar  spacings  of  the  S2  silicide.  Moreover,  all  of 
the  other  observed  interplanar  spacings,  with  the  excep¬ 
tion  of  0.222,  0.171,  0.146,  and  0.135  nm,  are  also  in 
good  agreement  with  those  cafculated  for  S2 .  These 
interplanar  spacings,  however,  match  well  with  those  of 
a-titanium.  Thus,  it  can  be  stated  that  S2  silicide  defi¬ 
nitely  exists  in  that  alloy. 

It  is  important  to  mention  here  that  Kotval  and  Calder"41 
also  have  characterized  silicides  in  the  complex  alloys 
Ti-5.\!  5Sn-2Zr-0.8Mo-0.7Si  and  Ti-4.5Al-2Sn-2Zr-3Ga- 
2Mo-0.5Si  as  S',  using  electron  diffraction  technique. 
Analysis  of  their  selected  area  diffraction  patterns  re¬ 
veals  that  the  silicide  in  alloy  Ti-4.5Al-2Sn-2Zr-3Ga-2Mo- 
0.5Si  could  be  indexed  either  as  S,  or  5'.  However,  the 
electron  diffraction  pattern  of  the  silicide  in  the  alloy 
Ti-5Al-5Sn-2Zr-0.8Mo-0.7Si  can  be  indexed  only  as  S' 
silicide. 

In  conclusion,  it  may  be  stated  that  all  of  the  inter¬ 
planar  spacings  for  the  silicide  in  alloy  Ti-6Al-3Sn-3Zr- 
0.5Si  reported  by  Antony151  are  accounted  for  in  terms 
of  the  typical  S2  silicide  and  the  matrix  a-titanium.  How¬ 
ever,  in  the  context  of  the  findings  of  Kotval  and  Calder,"41 
detailed  investigations  using  both  X-ray  as  well  as  elec¬ 
tron  diffraction  techniques  are  required  to  fully  resolve 
the  controversies. 
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Observation  of  a  Massive 
Transformation  from  a  to  y  in 
Quenched  Ti-48  At.  Pet  Al  Alloys 

PING  WANG,  G.B.  VISWANATHAN, 
and  VIJAY  K.  VASUDEVAN 

Owing  to  their  unique  combination  of  attractive  prop¬ 
erties,  such  as  lightness,  superior  strength,  creep  and  ox¬ 
idation  resistance,  alloys  based  on  the  intermetallic 
compound  TiAl  have  been  of  much  interest  in  recent  years 
as  structural  materials  for  elevated  temperature  aero¬ 
space  applications."1  The  alloys  developed  to  date  that 
show  the  most  promise  with  regard  to  room- temperature 
ductility  and  elevated  temperature  strength  are  based  on 
the  Ti-48A1  composition  (compositions  in  the  text  refer 
to  atomic  percent),  with  additions  of  elements  such  as 
V,  Cr,  Nb,  Mn,  etc.,  singly  or  in  combination.12'61  All 
of  these  alloys  are  two-phase,  with  a  microstructure 
composed  of  either  fully  lamellar  grains  of  Ti3Al  +  TiAl 
(a2  +  y)  or  duplex  primary  y  plus  lamellar  grains,  de¬ 
pending  on  processing  and  heat  treatment.141  The  prop¬ 
erties  of  these  alloys  are  quite  sensitive  to  microstructure 
and,  therefore,  to  heat  treatment;  duplex  structures  show 
higher  tensile  ductilities,  for  example,  up  to  2.2  pet  in 
binary  Ti-48A1|3-41  and  up  to  4  pet  in  ternary  or  quarter- 
nary  alloys, ,3_6!  whereas  fully  lamellar  structures  have 
poor  ductility.14"71 

The  lamellar  ( a2  +  y)  structure  was  first  studied  by 
Blackburn,181  who  reported  the  formation  of  alternate  la¬ 
mellae  of  y  and  a,  with  an  orientation  relationship  in¬ 
volving  {111}  y//{0001}  a2  and  (Oil)  y//{  1 1 20)  a2  and 
suggested  that  the  y  may  form  by  a  stacking  fault  mech¬ 
anism.  More  recently,  several  studies  of  phase  equilibria 
and  transformations  in  two-phase  Ti3Al  +  TiAl  alloys 
have  appeared.19"171  Of  particular  significance  are  the  re¬ 
sults  of  these  investigations  showing  decidedly  that  the 
(/ 3  +  y)  two-phase  field  shown  in  previous  Ti-Al  phase 
diagrams"81  does  not  exist  and,  instead,  the  (a  +  y)  phase 
field  extends  up  to  the  peritectic  temperature;  i.e.,  the  a 
region  exists  up  to  1450  °C  and  50  at.  pet  Al.  Recently, 
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several  investigators'1'*  have  studied  the  lamellar 
structure  and  reported  the  observation  of  different  types 
of  orientation  and  twin  variants  and  interfaces  of  TiAI. 
It  has  been  suggested'4 171  that  in  alloys  containing  up  to 
—42A1.  the  lamellar  structure  is  formed  by  the  reaction 
a  — 1 ►  as  — »  as  +  y,  whereas  at  higher  concentrations, 
the  y  lamellae  are  formed  directly  from  the  a  by  possibly 
a  stacking  fault-like  mechanism,  and  subsequently,  the 
a  orders  to  a2  to  produce  the  same  morphology  of  la¬ 
mellae  as  observed  at  lower  concentrations.  It  might  be 
expected  that  the  transformation  of  the  a  phase  would 
be  sensitive  to  cooling  rate.  Other  than  studies  on  the 
effect  of  rapid  solidification,1"1151  there  have  been  lim¬ 
ited  studies  of  the  effect  of  cooling  rate  on  the  transfor¬ 
mation  of  the  disordered  «  phase  in  high-Al  alloys.  The 
purpose  of  this  article  is  to  report  on  such  an  investi¬ 
gation  in  a  Ti-48AI  alloy.  In  particular,  the  observation 
of  a  massive  transformation  from  a  to  y  is  highlighted 
and  documented  for  the  first  time. 

A  Ti-48AI  alloy  (#1)  was  prepared  by  arc  melting 
commercial-purity  Ti  and  99.999  pet  pure  Al.  The  alloy 
button  was  homogenized  at  1150  °C  for  50  hours  fol¬ 
lowed  by  slow  cooling.  A  second  Ti-48A1  (#2)  cast  ingot 
was  obtained  from  The  Duriron  Company,  Inc.,  Dayton, 
OH,  hot  isostatically  pressed  at  1200  °C  for  2  hours, 
sealed  in  a  Ti-6AI-4V  can  under  vacuum,  and  extruded 
at  1360  °C  at  a  ratio  of  6:1.  The  Al  contents  after  ho¬ 
mogenization  in  the  two  alloys  were  determined  to  be 
47.20  and  47.86  at.  pet,  respectively.  Caution  was  taken 
to  minimize,  as  far  as  possible,  the  reaction  of  the  sam¬ 
ples  during  heat  treatment  with  the  atmosphere  and  the 
consequent  ingestion  of  oxygen  and  nitrogen.  Cut  sam¬ 
ples  were  placed  in  quartz  capsules,  evacuated  to  a  pres¬ 
sure  of  1(T6  torr  repeatedly  after  back-filling  with  argon, 
and  finally  sealed  under  argon.  In  addition,  Ti  foils,  placed 
separately  inside  the  capsules  well  away  from  the  sam¬ 
ples,  were  heated  to  about  800  °C  after  encapsulation  to 
getter  any  remnant  oxygen  and  nitrogen  in  the  back-filled 
argon.  Another  scheme  for  protection  used  consisted  of 
coating  the  samples  with  a  high-temperature  glass-ceramic 
(CORNING*  Glass  7913)  prior  to  direct  placement  in 

‘CORNING  is  a  trademark  of  Coming  Glass  Works.  Coming,  NY 


the  furnace.  Most  of  the  results  presented  herein  utilized 
the  latter  scheme.  Protected  samples  were  heated  to  the 
completely  a  region,  namely,  1400  °C,  for  various  times 
ranging  from  15  minutes  to  an  hour  and  then  cooled  at 
various  rates,  namely,  furnace  cooled  (FC),  air  cooled 
(AC),  oil  quenched  (OQ).  water  quenched  (WQ).  and 
iced  brine  quenched  (IBQ).  Examination  of  the  micro¬ 
structures  revealed  no  significant  effects  of  the  solution- 
izing  time  on  the  transformation  modes.  Finally,  to. address 
the  issue  of  the  possible  effect  of  oxygen  and  nitrogen, 
some  samples  were  solution  treated  in  air  for  various  times 
and  cooled  at  various  rates.  The  microstructures  of  these 
were  essentially  identical  to  those  that  received  protec¬ 
tion.  The  oxygen  and  nitrogen  contents  after  these  var¬ 
ious  treatments  were  analyzed  and  are  presented  later. 
Analysis  of  the  Al  contents  after  heat  treatment  indicated 
that  a  smali  loss,  typically  <0.25  at.  pet.  had  occurred. 

Optical  microscopy  observations  were  conducted  on 


samples  mounted,  polished,  and  etched  with  Kroll's  re¬ 
agent.  The  microstructures  were  also  observed  in  a  JEOL 
JXA-840  scanning  electron  microscope  (SEM)  utilizing 
back-scattered  electrons  for  obtaining  atomic  number 
contrast.  Fully  quantitative  electron  microprobe  analysis 
of  the  compositions  of  phases  was  conducted  in  the  same 
SEM.  Thin  foils  for  transmission  electron  microscopy 
were  prepared  from  the  heat-treated  samples  by  the  twin- 
jet  technique  using  a  solution  consisting  of  300  ml  meth¬ 
anol.  180  ml  butyl  alcohol,  20  ml  perchloric  acid,  and 
10  ml  hydrochloric  acid,  at  a  voltage  of  —70  V.  current 
of  20  to  30  mA,  and  a  temperature  of  -40  °C.  Obser¬ 
vations  of  the  foils  was  conducted  in  a  PHILIPS*  CM20 

'PHILIPS  is  a  trademark  of  Philips  Electronic  Instruments 
Corporation.  Mahwah.  NJ. 


transmission  electron  microscope  (TEM)  operated  at 
200  kV,  utilizing  bright-field  (BF),  dark-field  (DF),  and 
selected  area  diffraction  (SAD)  imaging  modes. 

The  optical  microstructurcs  of  samples  of  alloy  1  cooled 
at  different  rates  from  the  completely  a  region  are  shown 
in  the  micrographs  in  Figures  1(a)  through  (0-  Aside  from 
the  characteristic  fine  lamellar  structure,14 13171  another 
prominent  feature  present  in  the  FC  and  particularly  the 
AC  micrographs  (Figures  1(a)  and  (b))  is  a  structure  that 
appears  in  a  Widmanstatten  morphology,  in  three  to  four 
orientations  within  a  single  grain.  The  AC  sample  also 
shows  the  presence  of  deeply  etched,  mottled  regions. 
As  shown  later,  both  of  these  structures  were  established 
to  be  y.  A  dramatic  change  in  the  microstructure  is  ob¬ 
served  at  the  higher  cooling  rates.  A  remarkable  feature 
seen  in  the  micrographs  of  samples  that  were  quenched 
in  oil,  water,  and  iced  brine  (Figures  1(c)  through  (e)) 
is  the  presence  of  dark,  deeply  etched,  patchy  regions  in 
a  featureless  matrix.  These  regiot  s  appear  to  have  nu¬ 
cleated  largely  at  the  prior  a  grain  boundaries  and  then 
spread  rapidly  into  the  grains  in  a  spheroidal  and  massive 
fashion.  Often  these  dark  regions  were  also  found  within 
grains,  at  some  distance  from  the  grain  boundaries 
(Figure  1(e)).  The  interface  between  the  dark  and  light 
regions  varies  from  smooth  to  jagged  from  one  region 
to  another  and  appears  to  be  incoherent,  although  this 
has  not  been  established  positively.  The  volume  fraction 
of  the  dark  regions  did  not  appear  to  change  appreciably 
with  an  increase  in  cooling  rate. 

To  ascertain  whether  the  dark  regions  formed  only  upon 
quenching  from  the  fully  a  region,  another  sample  was 
heated  to  1400  °C,  slowly  cooled  to  1360  °C,  i.e..  into 
the  a  +  y  region.  and  WQ.  The  micrograph  correspond¬ 
ing  to  this  treatment,  shown  in  Figure  1(F).  appears  prac¬ 
tically  identical  to  those  from  samples  quenched  rapidly 
from  above  the  a  transus.  It  should  be  noted,  however, 
that  a  few  a  +  y  lamellae,  primarily  restricted  to  the 
grain  boundaries,  were  formed  during  the  initial  cooling 
from  1400  °C.  whereas  the  holding  at  1360  °C  prior  to 
the  quench  did  not  produce  any  primary  y  grains,  so  that 
the  composition  of  the  a,  for  the  most  part,  should  have 
remained  the  same  as  the  alloy  composition.  One  pos¬ 
sible  reason  for  the  absence  of  the  primary  y  grains  is 
related  to  the  fact  that  the  nucleation  of  these  is  de¬ 
pendent  on  the  prior  grain  size  and  holding  time  and 
both  the  large  grain  size  of  the  as-cast  +  homogenized 
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Fig.  1  —Optical  microstructures  of  Ti-47  2AI  alloy  held  at  1400  JC  for  1  h  and  cooled  at  various  rates  («)  1C.  I/O  AC.  (<•)  OQ  (</)  \vg 
(e)  IBO-  and  I  f )  slow  cooled  from  I4(X)  C  to  1760  ’C.  held  15  min.  followed  by  a  water  quench 


samples  (500  through  1000  /tm|  and  short  time  of  15 
minutes  diminished  the  opportunity  for  nudeation. 

In  order  to  further  understand  the  microstructures,  TEM 
observations  of  the  samples  cooled  at  the  different  rates 
were  conducted  The  light,  featureless  regions  in  the 
samples  quenched  at  high  cooling  rates  were  established 
to  be  <»..  as  shown  in  the  OF  micrograph  and  inset  dif¬ 
fraction  pattern  in  Figure  2  taken  from  the  WQ  sample: 
the  presence  of  antiphase  domains,  as  well  as,  on  oc¬ 
casion.  stacking  faults,  can  also  be  noted.  The  dark  re¬ 
gions  in  the  high  cooling  rate  samples  were  determined 
to  be  y.  An  example  is  shown  in  Figure  3  by  the  BF 


micrograph  and  SAD  patterns  corresponding  to  several 
y  /one  axes,  taken  from  the  WQ  sample.  The  lattice  pa¬ 
rameters  of  the  y  phase  were  determined  to  be  a  = 
4.HX)  A  and  c  =  4.138  A,  c/a  of  1009.  which  is 
slightly  more  cubic  than  the  c/a  of  1.016  for  stoichio¬ 
metric  y.  A  complex  subsiructure  consisting  of  antiphase 
boundaries  (APBs).  stacking  faults,  microtwins,  and  dis¬ 
locations  in  the  y  phase  can  be  seen  in  Figure  3(a).  Hlec- 
tron  diffraction  analysis  of  many  y  grains  showed  the 
occurrence  of  all  six  variants  Aside  from  face-centered 
cubic  (fee)  twin-related  variants,  numerous  contacting 
variants  in  perpendicular  twin  relation,  i.c..  mutually 
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Fig.  2  —  DF  TEM  micrograph  and  inset  1 12 1 .1 1  SAI)  pattern  showing 
the  microstructure  of  a,  in  the  Ti-47  2AI  alios  held  at  1400  for 
I  h  and  WQ 


perpendicular  c-axes,  arc  observed.  In  fact,  many  of  the 
boundaries  between  grains  were  observed  to  be  inter¬ 
variant  boundaries.  These  observations  suggest  that  a  di¬ 
rect  transformation  from  a  to  y  had  occurred,  perhaps 
in  a  massive  fashion.  Pools  of  the  lamellar  (o:  +  y)  grains 
were  also  observed  distributed  sparsely  between  the  y 
grains.  These  presumably  formed  due  to  recalescence 
effects. 

In  order  to  establish  that  indeed  a  massive  transfor¬ 
mation  had  occurred,  the  samples  were  observed  by  back- 
scattered  electron  (BSE)  imaging,  and  the  chemical 


Fig  3  — TEM  micrographs  of  Ti-47  2AI  alloy  held  at  1400  °C  for 
I  h  and  WQ  showing  massively  transformed  y  nil  BE  micrograph 
and  (b)  Ihrough  (el  are.  respectively.  1 101).  |II0|.  |00l|.  (010)  y 
SAD  patterns 


composition  of  the  phases  was  determined  by  micro¬ 
probe  analysis.  Thus.  BSE  micrographs  of  the  WQ  and 
IBQ  samples  are  shown  in  Figures  4(a)  and  (b),  respec¬ 
tively.  The  darker  and  lighter  contrast  regions  corre¬ 
spond  to  y  and  os,  respectively.  Although  the  observation 
of  contrast  differences  between  the  two  phases  indicates 
a  difference  in  chemical  composition,  a  microprobe  scan 
analysis  of  the  a:  and  y  regions  extending  through  the 
interface,  over  a  distance  of  ~7(X)  /am  at  —  l()-gm  in¬ 
tervals  (Figure  5).  reveals  that  the  compositions  of  the 
two  phases  are  practically  identical.  Another  feature  ob¬ 
servable  in  the  SEM  micrographs  (Figures  4(a)  through 
(b))  is  that  the  interface  seen  in  the  optical  micrographs 
(Figures  1(c)  through  ( e ) )  often  exhibits  ledge  features 
and  planar  and  crystallographic  facets,  implying  that  it 
may  actually  be  partially  coherent,  although  this  has  not 
been  established  positively. 


Fig  4  Back- scattered  SEM  micrographs  of  Ti -47  2AI  alloy  mi  WQ 
and  ( b)  IBQ  after  a  1400  ’C.  I  -h  treatment  The  dark  and  light  regions 
correspond  to  y  and  o-,  respectively 
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DISTANCE  FROM  START  OF  SCAN  (MICRONS) 

Hg  5  —  Micrnpmhc  scan  showi.ic  that  the  composition  tat  pet  Ah 
of  massive  y  and  «•  phases  in  the  Tt-47  2AJ  alloy  held  at  I4DO  C 
tor  \  h  and  WQ  are  practicallv  identical 

Several  points  with  regard  to  the  BSE  micrographs 
should  be  noted.  Both  the  microprobc  analysis  and  the 
BSE  imaging  were  conducted  on  polished,  unetched 
samples.  The  BSE  micrographs  in  Figures  4(a)  and  (b) 
appear  to  indicate  that  the  Al  content  of  the  y  phase  is 
slightly  higher  than  that  of  the  However,  the  micro¬ 
probe  analysis  (Figure  5)  does  not  show  any  such  con¬ 
sistent  trend  that  might  support  this.  It  was  practically 
impossible  to  distinguish  between  the  a:  and  massive  y 
regions  under  normal  contrast  settings  in  the  SEM.  At 
75  pet  of  the  maximum  setting,  differences  in  contrast 
became  just  discernible,  becoming  clearer  at  the  maxi¬ 
mum  setting.  It  is  possible  that  the  contrast  differences 
are  produced  by  a  channeling  effect  because  of  the  dif¬ 
ference  in  orientations  of  the  two  phases.12'1  Were  this 
the  case,  rotation  or  tilting  of  the  sample  should  have 
produced  changes/reversal  in  contrast,  which  was  not 
observed  Based  on  these  observations,  it  appears  that 
the  contrast  differences  seen  in  the  BSE  images  may  be 
associated  with  the  differences  in  crystal  structures  and. 
possibly,  density  and  defect  structures  of  the  two  phases, 
although  additional  work  is  required  to  determine  the  ex¬ 
act  cause.  In  any  case,  it  is  clear  that  the  compositions 
of  the  a:  and  y  phases,  within  the  experimental  accuracy 
of  the  electron  probe  microanalysis  technique,  are  nearly 
the  same,  which,  together  with  TEN1  results,  provides 
strong  evidence  for  a  massive  transformation  from  a  to 
y  at  high  cooling  rates 

To  ascertain  whether  the  transformations  were  sensi¬ 
tive  to  prior  grain  size,  the  extruded  Ti-4XAI  alloy  (#2) 
was  heated  to  1400  C  and  cooled  at  various  rates  Com¬ 
pared  with  the  unextnided  alloy  (#1)  with  a  grain  xi/e 
of  5<X)  to  l(XX)  /am.  the  as-extruded  microstructure  is 
duplex  (Figure  6(a)).  consisting  ol  primary  y  and  la 
mellar  la.  *•  y)  grains  with  grain  sizes  of  40  and 
-1(H)  /am,  respectively.  The  AC  microstructure  of  the 
latter  (Figure  Ob)  appears  with  a  larger  fraction  of  the 
Widmanstatten  structure  compared  w  ith  that  of  the  unex¬ 
truded  AC  sample  (Figure  1(b)).  The  WQ  sample 
(Figure  6(c))  appears  to  have  massively  transformed  to 
y  practically  throughout  the  complete  cross  section  of 
the  sample  A  feathery  morphology  for  the  massive 
product  was  also  observed 

A  BF  TEM  micrograph  of  the  AC  extruded  sample  is 
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shown  r.i  Figures  7ia>  The  large  feature  present  in  the 
center  .if  the  micrograph  corresponds  to  (he  Widmanstatten 
structure  seen  in  the  optical  micrographs  m  Futures  Kb) 
and  6ib)  By  means  ot  electron  diffraction  (Figure  7(b)). 
this  structure  was  determined  to  be  y  Csuallv.  three  or 
lour  variants  of  y  were  observed  to  be  present  within  a 
Widmanstatten  colony;  a  very  thin  layer  u|  rr,  could  also 
be  detected  The  micrograph  in  Figure  7(.u  ;,iso  shows 
the  intragrain  (o.  >  y)  lamellar  structure .  as  established 
by  the  composite  |  I  I0|v  *  |1<UK  tw  in  //  |  j  T 2()| ,  SAD 
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flL’  7  1 1M  micrographs  of  Ii  4?  S6A1  extruded  allov  AC'  after  a 

1400  C,  I  h  treatment  i*i»  10  micrograph  showing  Widmanstatten 
v  and  («  *  vi  lamellar  structure.  i/>i  |  !<•)  |y  •  ||0||\  i\% in  SAD 

pattern  ot  Wulmanstatten  structure,  and  i<  i  |  I  H)}v  *  JlOljv  twin  ' 
1 1  1  20  jo  SAD  pattern  from  if*  •  v*  structure 


pattern  in  Figure  7(c)  I  ransmission  electron  micro 
graphs  and  SAD  patterns  taken  from  the  extruded  WQ 
alloy  (Figure  X)  establish  that  the  n  has  practically  com¬ 
pletely  massively  transformed  to  v  In  terms  of  detect 
structures,  the  y  appears  with  features  very  similar  to 
those  observed  in  the  unextnided  WQ  alloy;  a  larger 
incidence  of  lee  twin-related  variants  was  noted  Hack- 
scattered  electron  images  showed  no  observable  differ¬ 
ences  in  contrast,  i  c  .  were  uniformly  dark  Based  on 
microprobe  scan  analysis,  the  composition  of  the  alloy 
was  determined  to  be  uniformly  near  the  value  ol 
47  86  at  pet  Al  determined  b\  wet  chemical  analvsis 
The  massive  transformation  has  been  extensively  stud 
led  in  the  past  '*  Prominent  examples  of  the  massive 
transformation  can  be  found  in  Cu  /n.  '  (ii  Al. 
Cu-Cia.1  '  Ag  \.  '  '  Fe-Ni.  and  I  i  -  A  1  alloys  l.s 

sentially.  a  massive  transformation  takes  place  by  dit- 
fusional  nuclealion  and  growth  aiul  involves  a  change  m 
crystal  structure  without  a  change  m  composition  : 
Generally,  the  massive  transiormation  is  observed  pro¬ 
vided  competing  equilibrium  and  other  inamelv. 
W  idmanstatten )  transformation  modes  are  suppressed 
l  sually  .  nuclealion  occurs  at  the  gram  boundaries  ol  the 
parent  phase,  and  growth  always  proceeds  by  a  ther¬ 
mally  activated  process  m  which  noncoopci.itnc  transfer 
ot  atoms  across  relatively  ineolieient  interfaces  oc 
curs  In  principle,  a  nijssice  transformation  can  be 

observed  ifl  any  system  whic  the  pertinent  phase  fields 
overlap  in  composition  Moreover,  it  is  possible  tor 

a  massive  reaction  to  be  considered  as  occnrmg  m  tile 
metastable  two  phase  field,  such  .is  that  below  a  mvari 
ant  reaction  temperature  itor  example,  euteetoul  tem- 
(X'rature)  Significant  penetration  into  an  equilibrium  or 


fhe  above  conditions  are  adeijuatelv  met  around  the 
two-phase  to  •  yi  region  in  the  Ti-AI  system  The  pres¬ 
ent  study  indicates  that  the  v  mav  form  from  the  a  in 
three  morphologies  depending  on  cooling  rate,  and 
therefore,  undercooling  At  low  value:,  the  lamellar 
morphology  prevails;  at  slightly  higher  rates,  the 
Widmanstatten  structure  appears,  and  .it  very  high  rates, 
the  massi.c  morphology  is  prevalent  I  he  electron 
microscopy  observations  and  the  observed  composition 
invariance  provide  strong  evidence  lor  the  massive  trans¬ 
iormation  from  n  to  y  The  extent  ol  the  Widmanstatten 
and  massive  morphologies  also  appears  to  de|Vtkl  on  pnor 
grain  si/e  For  instance,  compared  with  the  larger  grain 
si/e  unextruded  sample,  in  the  liner  grain  si/e  extruded 
sample,  the  volume  proportion  ol  both  morphologies  is 
observed  to  increase,  such  that  at  the  higher  cooling  rates, 
practically  the  entire  sample  is  transformed  to  massive 
y  One  explanation  ior  tne  difference  in  behavior  can  be 
forwarded  based  or.  a  grain  si/e  effect  Since  the  unex¬ 
truded  samples  of  alloy  1  are  characterized  bv  a  large 
grain  si/e  iXiH)  io  liHKi  /ami.  lewer  sites  tor  nueleation 
ol  the  massive  y  exist  Once  nucleated,  these  grow  at 
high  rales  until  the  temperutute  drop,  sufficiently  below 
I  I  25  (  to  encounter  the  extended  <«/.»  phase  boundaiv  . 
whereupon  the  transformation  ot  the  remaining  o  r.  o 
intervenes  and  arrests  further  progression  of  the  massive 
reaction  Another  possible  contribution  comes  limn  (he 
higher  M  content  ot  the  extruded  material,  which  might 
permit  the  massive  leaetion  io  piogrcss  tor  longer  limes 
before  the  o  ’««  boundaiv  is  encountered  I  he  relative 
c onli ibutioiis  ol  these  two  etlects  is  being  investigated 
Finally,  while  morphologically  the  Widmanstatten  and 
massive  y  structures  appear  different,  there  is  uncei 
laintv  in  distinguishing  between  the  two  in  terms  ol 
mechanisms  and  tcmpcialuic  dependence.  since  the 
composition  ol  the  Widmanstatten  >  has  not  been 
ascertained 

An  import. ini  question  iciates  to  the  temperaiure  ot 
mic  leution  of  ifie  massive  v  It  is  geneiallv  believed  that 
the  /  line,  i  <  .  the  loiiis  ,q  points  ,u  which  ;he  lice 
enerev  ol  llie  two  phases  o  equal.  lepie'enls  the  limit 
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for  the  massive  reaction  and  the  distance  between  the  Ta 
line  and  the  temperature  of  the  start  of  the  reaction  rep¬ 
resents  the  undercooling  and.  therefore,  the  driving  force 
for  the  massive  transformation.124  251  Hillert,132'  however, 
has  argued  that  the  massive  transformation  requires 
undercooling  below  the  solvus  line.  Another  point  con¬ 
cerns  the  nature  of  the  interface  between  the  parent  and 
the  massive  product.  In  this  context,  it  is  noteworthy  that 
evidence  for  both  incoherent  and  coherent  growth  during 
the  massive  transformation  is  a  subject  of  much  de¬ 
bate,1331  since  the  nature  of  the  massive /parent  interface 
varies  from  incoherent  to  semicoherent  from  system  to 
system.  Evidence  for  both  the  absence  and  presence  of 
definite  crystallographic  relationships  between  the  mas¬ 
sive  product  and  parent  exists.133  341  Even  so,  it  has  been 
argued  that  massive  transformations  must  always  involve 
special  crystallographic  relations  (between  the  parent  and 
the  product),  leading  to  coherent  nucleation  followed  by 
growth,  in  which  the  interface  changes  from  coherent  to 
semicoherent  and  moves  by  a  ledge  mechanism.1341  Our 
current  experimental  evidence  suggests  that  the  interface 
between  the  a2  and  massive  y  may  be  semicoherent. 
However,  determination  of  the  exact  nature  of  the  inter¬ 
face  by  the  TEM  will  require  foils  in  which  the  massive 
product,  the  parent,  and  the  separating  interface  are  pres¬ 
ent  simultaneously,  which  is  the  subject  of  ongoing  re¬ 
search.  All  of  these  aspects,  together  with  the  temperature 
and  composition  dependence  and  the  nucleation  and 
growth  mechanisms  of  the  different  morphologies  of  y, 
are  the  subject  of  ongoing  research. 

A  final  point  concerns  the  possible  role  of  oxygen  and 
nitrogen.  Various  studies  of  the  possible  effects  of  inter¬ 
stitials  were  carried  out  in  which  the  solution  treatment 
time  was  varied  from  1 5  minutes  to  an  hour  in  both  pro¬ 
tected  as  well  as  unprotected  samples.  The  interstitial 
contents  analyzed  following  some  of  these  treatments  are 


shown  in  Table  I,  from  which  it  is  clear  that  pickup  of 
interstitials,  particularly  oxygen,  is  a  function  of  time 
solutionizing  and  the  protection  scheme.  At  the  longer 
times  (1  hour),  a  fair  amount  of  oxygen  (—0.3  wt  pet) 
appears  to  have  been  ingested  into  the  samples  compared 
with  those  at  shorter  times  (~0.096  to  0.14  wt  pet).  (It 
is  thought  that  the  increased  oxygen  in  many  of  these 
samples  was  caused  by  reaction  with  the  glass  ceramic 
coating  or  by  contact  with  the  quartz  tubes  used  to  pro¬ 
tect  the  samples.)  Even  so,  the  microstructures  of  these 
samples  cooled  at  various  rates,  including  those  cooled 
at  high  rates,  with  regard  to  the  mode  and  morphology 
of  the  transformation  products,  remained  essentially  the 
same  (Table  I),  suggesting  that  the  interstitial  contents 
do  not  exert  a  dominant  influence  on  the  transformation. 
Lastly,  an  important  point  to  also  be  noted  in  Table  I  is 
that  chemical  analysis  of  the  interior  of  the  samples  after 
removing  ~2  mm  of  the  outer  surface  revealed  that  the 
pickup  of  oxygen  was  confined  to  the  outer  surface, 
whereas  the  interiors  were  practically  unaffected. 

In  view  of  the  foregoing,  it  is,  nevertheless,  interest¬ 
ing  to  examine  the  potential  role  that  oxygen  could  play. 
In  addition  to  the  48A1  alloy,  studies  were  also  con¬ 
ducted  on  a  nominal  Ti-50A)  alloy  cooled  at  various  rates 
from  the  a-region.  The  behavior  of  the  latter  was  similar 
to  the  48A1  alloy;  both  the  Widmanstatten  and  massive 
y  morphologies  were  observed  depending  on  cooling  rate. 
Very  recently,  two  publications,  the  first  by  Jones  and 
Kaufmann1171  on  phase  equilibria  and  microstructures  in 
TiAl-base  alloys  and  the  second  on  the  effect  of  cooling 
rate  on  transformations  of  the  a  phase  in  Ti-48A1- 
2Mn-2Nb  alloys  by  McQuay  et  a/.,1351  have  appeared. 
Of  particular  interest  to  the  present  work  are  the  results 
in  the  first  on  a  50A1  alloy  solution  treated  at  1400  °C 
for  2  hours  and  1450  °C,  for  1  hour  and  then  WQ. 
The  microstructures  following  these  treatments  are  very 


Table  I.  Oxygen,  Nitrogen  Contents,  and  Microstructure  as  a  Function  of  Heat  Treatment  of  Ti-47.86A1  Alloy 


Heat  Treatment* 

Wt  Pet  Oxygen 

Wt  Pet  Nitrogen 

Microstructure** 

As-cast-hipped 

0.049  . 

0.006 

L  +  TV 

H  +  1400  °C.  20  min,  FC 

0.096 

0.009 

L  +  little  y. 

H  +  1400  °C,  1  h.  FC 

0.195 

0.009 

L  +  little  yH 

H  +  1400  °C.  20  min,  AC 

0.104 

0.008 

L  +  y. 

H  +  1400  °C,  1  h.  AC 

0.255 

0.011 

L  +  y. 

H  +  1400  °C,  1  h,  AC’ 

0.040 

0.005 

L  +  y„ 

H  +  1400  °C,  20  min,  OQ 

0.133 

0.008 

L  +  y.  +  y„ 

H  +  1400  °C.  1  h.  OQ 

0.311 

0.008 

L  +  y.  +  ym 

H  +  1400  °C.  20  min,  WQ 

0.119 

0.009 

mostly  y„  +  little  a2  +  L 

H  +  1400  °C  1  h.  WQ 

0.285 

0.008 

mostly  y„  +  little  a,  +  L 

H  +  1400  °C,  1  h,  WQ’ 

0.036 

0.004 

mostly  y„  +  little  a;  +  L 

As-extruded  (X) 

0.065 

0.033 

L  +  yp 

X  +  1400  °C.  20  min,  AC 

0.140 

0.033 

L  +  y„ 

X  +  1400  °C,  1  h,  AC 

0.319 

0.158 

L  *-  y„ 

X  +  1400  °C.  20  min,  WQ 

0.112 

0.018 

mostly  y„  +  little  a: 

X  +  1400  °C.  1  h,  WQ 

0.270 

0.119 

mostly  ym  +  little  a. 

H  +  1400  °C.  20  min  (air),  WQ 

0.297 

0  029 

Q;  +  y„ 

H  +  1400  °C,  1  h  (air),  WQ 

0.900 

0.180 

as  +  ym 

H  +  1400  °C,  1  h  (air),  WQf 

0.052 

0  005 

as  +  y„ 

*H  =  cast  +  hipped  (1200  °C,  2  h.  105  MPa);  X 

=  as-extruded;  FC  = 

furnace  cooled;  AC  =  air 

cooled;  OQ  =  oil  quenched;  and  WQ  = 

water  quenched 

**L  =  lamellar  (a;  +  y);  yp.  y..  and  y„  are,  respectively,  primary.  Widmanstatten,  and  massive  y 

Analysis  of  intenor  of  samples  after  removing  — 

2  mm  of  outer  surface 
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similar  to  those  observed  by  us  in  the  48  and  50A1  al¬ 
loys,  both  the  Widmanstatten  and  massive  morphologies 
could  be  noted.  Indeed,  it  is  significant  that  even  though 
the  measured  oxygen  and  nitrogen  contents  after  heat 
treatment  in  their  study  of  the  50A1  alloy  differ  appre¬ 
ciably  for  the  two  temperatures  (oxygen  0.051  and 
0.163  wt  pet,  nitrogen  0.021  and  0.09  wt  pet,  respec¬ 
tively).  the  microstnictures  are  essentially  the  same.  The 
microstructures  observed  at  low  (FC,  AC)  and  high  cool¬ 
ing  rates  (OQ,  WQ)  by  McQuay  et  al.  in  alloys  with  an 
oxygen  content  of  0.08  wt  pet  after  heat  treatment  are 
also  very  similar.  Indeed,  the  occurrence  of  massive  y 
at  high  cooling  rates  in  these  alloys  has  been  confirmed.1361 

Finally,  in  our  study  of  the  unextruded  samples  so¬ 
lution  treated  in  air  at  1400  °C  for  times  from  15  minutes 
to  an  hour,  the  massive  reaction  was  observed  to  be 
slightly  lesser  in  extent  near  the  outer  surface  compared 
to  those  in  the  protected  samples  (Table  I),  This  might 
suggest  that  the  high  oxygen  and  nitrogen  ingested  into 
the  surface  regions  of  the  former  act  to  inhibit  rather 
than  permit  the  massive  reaction  from  a  to  y.  Alter¬ 
nately,  it  is  possible  that  sufficient  loss  of  Al  occurred 
to  prevent  the  massive  transformation,  since  alloys  with 
Al  contents  less  than  —45  at.  pet  WQ  from  the  fully 
a-phase  region  transformed  completely  to  a2  and  showed 
no  evidence  of  massive  y.  However,  analysis  of  the  Al 
content  after  heat  treatment  in  both  the  surface  and  in¬ 
terior  of  the  samples  showed  no  appreciable  change  from 
the  value  before  heat  treatment.  If  indeed  the  ingestion 
of  interstitial  oxygen  is  able  to  produce  such  dramatic 
changes  in  the  transformation  mode,  that,  in  itself,  poses 
some  rather  intriguing  questions  with  regard  to  how  such 
a  change  is  brought  about,  since  oxygen  is  known  to 
stabilize  the  a  phase.  There  is  one  instance,  namely,  pure 
Ti.  where  because  of  the  addition  of  oxygen,  a  massive 
transformation  from  (i  to  a  has  been  postulated  to  occur 
prior  to  the  martensitic  transformation.1371  Thus,  one  po¬ 
tential  way  in  which  oxygen  could  have  an  effect  is  by 
altering  the  a  and  y  phase  boundaries  such  that  the  over¬ 
lap  in  compositions  of  the  two  phases  at  different  tem¬ 
peratures  is  increased,  thereby  permitting  the  massive 
reaction  to  occur  over  a  wider  range  of  Al  contents.  Al¬ 
though  the  possibility  that  even  very  low  levels  of  ox¬ 
ygen  can  produce  such  changes  cannot  be  overruled,  our 
current  experimental  evidence,  together  with  the  results 
of  others117-35'  discussed  in  the  foregoing,  is  consistent 
with  the  suggestion  that  the  massive  transformation  is  a 
real  effect,  not  some  chance  consequence  of  ingestion  of 
interstitials,  particularly  oxygen.  Additional  work  on  low 
oxygen  alloys  solutionized  and  quenched  in  vacuum  has 
essentially  confirmed  this. 
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4.0  SUMMARY 


A  series  of  hot -compress ion  tests  were  conducted  on  coarse-grained 

wrought  alloys.  Monel  K-500  and  alloy  718.  The  most  interesting  results  came 

from  tests  conducted  with  a  deformation  level  of  30X,  a  temperature  of  1125°C, 

and  a  strain  rate  of  1.0  s'1.  A  major  difference  in  recrystallization 

behavior  between  the  two  alloys  was  observed.  In  particular,  the  overall 

recrystallization  rate  for  Monel  K-500  was  higher  than  that  for  alloy  718; 

however,  the  nucleation  rate  at  boundaries  in  the  Monel  alloy  was  extremely 

variable,  and  this  feature  gave  rise  to  large  "patches"  In  the  microstructure 

which  remained  unrecrystallized  even  after  about  50X  of  the  microstructure  had 

recrystallized.  Some  "patches"  were  larger  than  1.5  millimeters,  which  would 

correspond  to  a  size  greater  than  100  millimeters  In  typical  castings  of  this 

alloy.  It  Is  postulated  these  phenomena  are  responsible  for  the  observation 

of  very  coarse  grains  in  billets  of  Monel  K-500.  In  contrast,  alloy  718 

exhibits  a  well-defined  "necklace"  structure  which  results  from  a  rather 

* 

uniform  nucleation  rate  at  almost  all  grain  boundaries.  As  a  result,  the 
microstructures  of  billets  of  alloy  718  should  be  uniformly  fine  grained. 
Recommendations  for  improving  the  uniformity  of  microstructure  in  billets  of 
Monel  K-500  are  included  in  the  body  of  this  report. 
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ABSTRACT 


The  primary  use  of  in  situ  seismic  measurements  in  soil  is  to  evaluate 
wave  velocities  and,  hence,  to  characterize  elastic  moduli.  However, 
material  damping  characteristics  of  soil  are  also  important  but  are  rarely 
measured  in  situ.  The  crosshole  seismic  method  is  shown  to  be  a  very 
useful  method  for  performing  attenuation  (damping)  measurements  at  small 
strains. 

To  initiate  developing  a  field  method  of  measuring  in  situ  attenuation 
characteristics,  testing  equipment  and  a  test  site  developed  and 
characterized.  Some  of  important  variables  associated  field  implementation 
of  damping  measurements  were  investigated.  Important  variables  are:  1 . 
the  frequency  content  of  the  waveform,  2.  the  ratio  of  the  wavelength  to  the 
source-receiver  spacing,  3.  the  frequency  characteristics  of  receivers  in 
boreholes  and  4.  the  quality  of  casing  and  grouting.  The  use  of  windowing 
time-domain  records  to  minimize  interference  from  reflected  and  refracted 
waves  and  the  calculation  of  spectral  ratios  in  the  frequency  domain  are 
necessary  approaches  to  data  processing  in  these  tests.  Also,  difference  in 
frequency  characteristics  between  two  receivers  in  boreholes  play  a  major 
role  in  damping  measurement. 

•» 

Recommendations  are  suggested  for  further  study  as  follows: 

1 .  the  test  site  needs  to  be  fully  characterized  by  using  other  in  situ 
testing  technique  such  as  Standard  Penetration  Test,  Cone 
penetration  Test  and  Pressuremeter  and  can,  thus  be  designated  as 
a  one  of  national  geotechnical  experimentation  sites 

2.  The  effect  of  casing  and  grouting  on  damping  measurements 
should  be  investigated.  That  .can  be  done  by  drilling  holes  and 
burying  receivers  at  specific  depths  without  casing  and  grouting. 

3.  The  Utilization  of  SH-waves  on  damping  measurements  is 
recommended.  Further  refinement  of  the  in-hole  source  with  a 
solenoid  developed  in  this  project  is  expected  to  make  it  possible. 
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Executive  Summary 

The  following  semiannual  progress  report  is  submitted  in  partial  fulfillment  of  the 
requirements  of  the  Air  Force  Engineering  Foundation  Grant  RI-B-90-06.  The  final  report 
is  due  October  15, 1991. 

Introduction 

The  specific  aim  of  this  research  is  to  develop  an  automated  technique  for 
quantifying  concrete  microcracking  using  the  principles  of  quantitative  stereology 
(Underwood,  1970).  The  goal  is  to  investigate  the  automation  of  the  oriented  secant 
method  either  by  a  digitized  system  or  by  using  an  image  analyzer,  and  to  determine  the 
effectiveness  of  this  computerized  procedure  in  characterizing  concrete  microcracking. 

Data  on  concrete  cylinder  specimens  subjected  to  uniaxial  compression  obtained  by 
Cornell  University  investigators  in  the  mid-1960’s  is  being  used  to  establish  the  validity  of 
this  automated  method.  (Hsu,  et  aL,  1963).  The  Cornell  investigators  used  a  cumbersome 
and  difficult  procedure  to  measure  and  relate  concrete  microcracking  to  uniaxial  stress- 
strain  curves.  The  fully  computerized,  automated  method  being  developed  uses  the 
oriented  secant  method  to  quantify  internal  microcracking.  This  automated  procedure 
will  enable  investigators  to  efficiently  characterize  concrete  microcracking  with  respect  to 
density  and  orientation.  For  more  information  concerning  the  Cornell  work  and  the 
oriented  secant  method  please  see  original  proposal  entitled  "Quantitative  Stereology  of 
Concrete  Microcracking." 
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COMPRESSIVE  FAILURE  DUE  TO  KINKING  OF  FIBROUS  COMPOSITES 

by 

Dimitri*  C.  Lagoudas  and  Ahmed  M.  Saleh 
Department  of  Civil  and  Environmental  Engineering 
and  Institute  Center  for  Composite  Materials 
Rensselaer  Polytechnic  Institute,  Troy,  NY  12180 


ABSTRACT 

This  work  addresses  the  important  engineering  problem  of  the  compressive  strength 
of  fibrous  laminated  composites,  bridging  the  micromechanics  relevant  to  compressive 
failure  due  to  fiber  kinking  with  the  global  characteristics  of  a  composite  structure. 
Simple  modeling  incorporates  the  influence  of  boundary  conditions,  laminate 
thickness  and  lay-up  configuration,  as  well  as  the  importance  of  the  micro-geometry. 
The  mechanism  of  failure  under  compression  is  assumed  to  be  micro  buckling  of  fibers 
that  localizes  and  leads  to  the  formation  and,  under  critical  conditions,  propagation  of 
kink  bands.  To  analyze  the  mechanics  of  kinking  and  calculate  the  critical 
compressive  stress,  a  steady  state  ldnk  propagation  model  is  utilized  for  the  layers 
under  compression  in  the  fiber  direction.  The  analysis  results  are  used  to  gain  a 
better  understanding  of  the  influence  of  the  macro-geometry  in  addition  to  the 
micro-geometry  and  the  fiber  and  matrix  material  properties.  The  predictions  of  the 
model  are  compared  with  experimental  results  for  carbon/thermoplastic 
unidirectional  laminates  under  direct  compression  and  cross-ply  laminates  under 
four-point  bending. 
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Friction  of  Granular  Material  Along  an  Interface: 
Mechanical  Modelling 


By  Samuel  G.  Paikowsky1,  Member ,  ASCE 


Abstract 

The  need  to  predict  the  friction  along  an  interface  between  a  solid  body  and 
granular  material  exists  in  many  different  disciplines.  However,  there  is  currently 
no  model  which  describes  or  explains  this  friction/shear  mechanism.  The 
controlling  parameters  of  this  mechanism  are  studied  through  a  mechanical  model 
enabling  to  predict  the  friction  coefficient  at  the  interface.  This  mechanical  model 
(S.G.  -  Sphere  in  the  Groove)  considers  the  micro  behavior  of  granular  material 
examining  a  single  rigid  grain  (sphere  or  ellipsoid)  supported  by  adjacent  particles 
(frictional  groove)  and  subjected  at  an  additional  contact  point,  to  normal  and 
traction  forces.  These  forces  can  result  from  a  moving  frictional  plane  (modeling 
the  body  in  contact  along  an  interface)  or  from  the  motion  of  an  adjacent  particle 
(modeling  shear  conditions  within  the  granular  material  itself). 

Based  on  observations  of  granular  $oil  behavior,  the  S.G.  model  assumes  that 
the  grain  will  roll  or  slide  along  the  groove  between  the  supporting  particles.  A 
limiting  equilibrium  analysis  is  applied  to  the  interface  and  interparticle  contacts, 
examining  the  possibility  of  sliding  or  rolling  in  both  locations.  A  companion  paper 
examines  the  S.G.  model  in  light  of  experimental  work  measuring  the  friction  of 
sands  along  interfaces.  These  analyses  are  then  followed  by  simplified  relations 
describing  the  coefficient  of  friction  between  quartzic  sand  and  solid  bodies. 


5  Ass't  Professor,  Dept,  of  Civil  Eng.,  University  of  Lowell,  Lowell,  MA  01854. 
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Conclusions  and  Suggestions  for  Future  Work 

An  instrument  has  been  developed  which  can  be  used  to  measure  spatial  concentration  in  an  air- 
water  mixture.  The  instrument  not  only  provides  the  capability  for  spatial  concentration  determination, 
but  also  allows  for  the  measurement  of  water  film  thickness  as  well  as  potential  for  flow  pattern 
identification  including  identification  of  the  transition  regions  between  flow  patterns.  The  instrument 
(CAC)  effectively  measures  the  spatial  concentration  and  film  thickness  of  water  and  air  in  the  active 
space  of  the  capacitive  sensor  in  nonobstructive  manner. 

Visually  observed  flow  regimes  (CCD  system  and  flow  visualization  system)  were  compared  with 
fluctuation  of  spatial  concentration,  where  the  frequency  distribution  (FDC)  and  die  cumulative  frequency 
distribution  (CFDC)  curves  of  spatial  concentration  were  the  characteristics  chosen  to  be  discriminators 
of  the  flow  pattern. 
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The  project  also  resulted  in  a  number  of  conclusions  relevant  to  instrument  development  and 

implementation: 

1.  Increasing  the  number  of  capacitive  sensors  may  provide  several  benefits  including  increased 
resolution  and  more  accurate  measurements,  and  identification  of  the  interpretative  algorithm  as 
well  as  control  of  a  calibration  status. 

2.  The  electrodes  should  be  approximately  in  the  range  of  0.1  to  2  diameters  in  length;  longer 
electrodes  will  yield  larger  averaging  in  the  output  while  shorter  lengths  will  result  in  some 
decrease  of  the  measured  delta  value.  For  average  cross-sectional  spatial  concentration  analysis, 
rapid  response  time  is  not  crucial.  However,  for  time  dependent  analysis,  as  in  the  measure  of 
spatial  concentration  fluctuations  (AC  component),  a  minimum  electrode  length  is  desirable  to 
minimize  spatial  averaging  in  this  parameter  value. 

3.  In  order  to  obtain  an  accurate  concentration  value  in  water-air  mixture,  the  correction  factor  for 
each  flow  pattern  should  be  carefully  determined  and  evaluated.  It  is  suggested  that  the 
mathematical  model  with  on-line  flow  pattern  detector  should  be  further  developed  in  order  to 
increase  accuracy. 

4.  Two  types  of  capacitive  sensors  have  been  employed  in  the  present  study  (types  I  and  III),  and 
the  sensors  seems  to  be  reliable  and  stable  enough  as  a  base  for  further  optimization  process  of 
the  sensors  for  the  measurement  of  spatial  concentration  and  film  thickness  of  a  liquid  in  a 
channel. 

5.  The  effects  of  spatial  concentration  and  the  temporal  and  spatial  distribution,  as  well  as  a 
configuration  of  a  capacitive  sensor  on  the  sensor  characteristic  are  significant. 

It  should.be  noted  tf  it  the  general  conclusions  drawn  regarding  the  capacitive  type  concentration  meter 

for  water-air  mixture  are  of  a  preliminary  and  a  exploratory  nature.  No  extensive  study  of  such 

parameters  as:  temperature,  conductivity,  size  and  configuration  of  electrodes,  dielectric  constant,  sensor 
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configurations,  etc...,  was  undertaken.  It  is  recommended  that  further  research  needs  to  be  made  as 
foliows: 

*  The  computed  spatial  concentration  is  somewhat  flow  regime  dependent.  Correction 
factors  could  be  incorporated  to  compensate  this  process. 

*  The  shape  of  the  F DC  possesses  a  high  potential  for  use  in  objective  determination  of 
flow  patterns  in  two-phase  flow.  It  should  be  noted  that  further  study  of  the  flow  pattern 
identification  would  seem  to  be  recommended. 

*  The  experimental  evidence  obtained  in  this  study  indicates  that  many  variables  effect  the 
capacitance  of  the  sensor  and  that  establishment  of  the  desired  value  may  be  difficult. 
It  is  suggested  that  calibration  test  need  to  be  developed  for  a  high  controllable  conditions 
for  such  variables  as  spatial  concentration  and  its  temporal  and  spatial  distribution,  and 
sensor  configuration. 
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Chapter  1.  Overview 


This  report  summarizes  the  results  of  the  research  program,  "Effects  of 
Reinforcement  Geometry  on  Inelastic  Deformation  in  Metal-Matrix  Composites,"  sponsored 
by  the  Air  Force  /  Engineering  Foundation  Grant  RI-B-90-10.  This  research  program 
consists  of  both  analytical  and  experimental  phases. 

The  analytical  phase  is  the  development  of  an  appropriate  micromechanics  model  to 
predict  the  elastic-plastic  response  of  metal-matrix  composites  with  ellipsoidal  inclusions. 
The  elastic-plastic  response  of  metal-matrix  composites  has  received  considerable  attention 
in  recent  years.  Due  the  complex  microstructure  of  these  materials,  analytical  modeling 
efforts  are  forced  to  idealize  the  geometry  in  order  to  define  solvable  boundary  value 
problems.  The  study  of  the  elastic-plastic  response  parallels  earlier  studies  of  the  elastic 
response.  Typical  approaches  for  the  elastic  properties  include  the  use  of  dilute 
approximations  (Dewey,  1947),  self-consistent  approximations  (Hill,  1965  and  Budiansky, 
1965),  idealized  geometries  (such  as  Hashin’s  composite  sphere  model  (1962),  Hashin  and 
Rosen’s  composite  cylinder  model  (1964),  and  Christensen  and  Lo’s  (1979)  three-phase 
model),  and  approximate  methods  based  on  the  so-called  Mori-Tanaka  (1973)  method. 
Many  of  these  models  make  use  of  Eshelby’s  (1957)  elasticity  solution  for  a  single  ellipsoidal 
inclusion  in  an  infinite  matrix.  In  the  present  study,  the  dilute,  self-consistent  and  three- 
phase  models  are  extended  to  study  the  composite  elastic-plastic  response.  Since  closed 
form  plasticity  solutions  (analogous  to  Eshelby’s  solution)  are  not  available,  numerical 
elastic-plastic  solutions  are  obtained  to  provide'  dilute  and  three-phase  estimates  of  the 
effects  of  reinforcement  geometry  on  the  composite  elastic-plastic  response.  A  generalized 
formulation  has  been  developed  to  allow  direct  comparisons  between  the  dilute,  self- 
consistent  and  three-phase  model.  For  the  case  of  spherical  inclusions,  a  Rayleigh-Ritz 
procedure  is  used  to  solve  the  boundary  value  problems  defined  by  each  model.  The  results 
of  the  three  models  are  compared  and  discussed  in  detail  in  Chapter  2  of  this  report  as  well 
as  in  the  recent  paper  by  Taggart  (1991).  It  is  concluded  that  the  three-phase  approach 
provides  reliable,  reasonable  estimates  of  the  macroscopic  response  as  well  as  details  of  the 
microstructural  deformation  mechanisms. 
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Since  the  dilute  approximation  of  the  elastic-plastic  response  is  a  straightforward  (yet 
numerically  intensive)  calculation,  the  effects  of  reinforcement  geometry  are  first  studied 
using  the  dilute  approach.  The  boundary  value  problem  defined  by  the  dilute  model  is  that 
of  a  single  inclusion  embedded  in  an  infinite  elastic-plastic  matrix  subjected  to  remote 
uniaxial  stressing.  To  solve  this  problem,  a  elastic-plastic  finite  element  analysis  program  was 
developed  and  applied  to  the  case  of  ellipsoidal  inclusions  of  various  aspect  ratios.  Note 
that  the  numerical  approach  allows  for  future  studies  of  arbitrary  axisymmetric  inclusion 
geometries.  While  the  results  of  the  dilute  model  are  expected  to  be  valid  only  for  very  low 
inclusion  volume  fractions,  it  is  anticipated  that  the  qualitative  trends  will  be  consistent  with 

improved  models.  Chapter  3  summarizes  the  results  of  the  development  and  application  of 

/' 

a  dilute  model  for  the  effects  of  inclusion  geometry  on  the  composite  elastic-plastic  response. 
These  results  are  also  presented  in  the  recent  paper  by  Adley  and  Taggart  (1991). 

In  chapter  4,  the  three-phase  model  studied  in  chapter  2  for  the  case  of  spherical 
inclusions  is  extended  to  study  the  effects  of  reinforcement  geometry  on  the  composite 
elastic-plastic  response.  For  the  case  of  ellipsoidal  inclusions,  the  application  of  the  three- 
phase  model  requires  additional  assumptions  which  are  not  required  for  the  case  of  spherical 
inclusions.  The  assumptions  are  discussed  and  evaluated  in  detail  in  chapter  4.  As  in  the 
case  of  spherical  inclusions,  the  dilute  and  three-phase  approximations  provide  similar  results 
at  low  inclusion  volume  fractions.  At  increased  volume  fractions,  the  three-phase  model  is 
shown  to  provide  a  stiffer  composite  response. 

The  experimental  phase  of  the  program  is  summarized  in  Chapter  5.  This  study 
consists  of  both  microstructural  and  mechanical  characterizations  of  commercial  extruded 
metal-matrix  composites  with  both  particulate  and  whisker  reinforcement  geometries.  The 
microstructural  characterization  consists  of  both  optical  and  scanning  electron  microscopy 
of  the  as-received  composite  microstructure.  Quantitative  analysis  of  these  micrographs 
provide  information  as  to  the  reinforcement  geometry,  aspect  ratio,  volume  fraction, 
orientation  and  distribution.  Additional  microstructural  characterization  consist  of  the 
examination  fracture  surfaces  and  internal  damage  associated  with  mechanical  loading. 
These  results  give  insight  into  the  critical  failure  mechanisms.  The  mechanical  tests  consists 
of  uniaxial  tension  tests  of  specimens  sectioned  both  parallel  and  normal  to  the  extrusion 
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direction  in  both  the  as-received  state  and  after  the  heat  treatment  recommended  by  the 
supplier.  The  elastic-plastic  response  for  these  composites  are  compared  to  each  other  and 
correlated  with  the  analytical  predictions.  It  is  concluded  that  for  the  case  of  particulate 
composites,  for  the  material  system  and  inclusion  volume  fraction  evaluated,  the  dilute  and 
three-phase  models  give  nearly  identical  estimates  of  the  composite  response.  The  predicted 
elastic  modulus  and  post-yield  strain  hardening  agrees  with  the  observed  behavior  but  the 
predicted  yield  stress  underestimates  the  experimental  yield  stress.  For  the  case  of  whisker 
reinforced  composites,  both  models  give  good  estimates  of  the  elastic  modulus  and 
overestimate  the  composite  yield  stress.  The  three-phase  model  gives  a  good  estimate  of 

the  post-yield  strain  hardening  behavior  while  the  dilute  model  underestimates  the  strain 

/ 

hardening.  The  discrepancies  between  the  experimental  and  theoretical  results  are 
attributed  to  the  models’  failure  to  consider  the  effects  of  processing  induced  residual 
stresses  and  internal  damage  mechanisms.  Future  studies  into  the  effect  of  residual  stresses 
and  damage  mechanisms  on  the  composite  macroscopic  response  are  planned. 


43 


RI-B-90-11 


Wake -Body  Interactions  and  the  Formation  of  Vortex  Shedding 

Behind  Bluff  Bodies 


Professor  Paul  J.  Strykowski 
Department  of  Mechanical  Engineering 
University  of  Minnesota 


44 


Wake-Body  Interactions  and  the  Formation  of  Vortex 
Shedding  Behind  Bluff  Bodies 


Abstract 

The  present  study  investigates  the  flow  around  a  circular 
cylinder  for  the  range  of  Reynolds  numbers  of  49  <  Re  <  300.  A 
thorough  investigation  was  done  on  the  boundary  conditions  such 
that  a  ‘universal’  relationship  between  the  Strouhal  number  and 
Reynolds  number  was  established.  The  boundary  conditions 
considered  included  flow  uniformity,  cylinder  vibration,  and 
cylinder  end  conditions.  The  effect  of  the  hot  wire  probe  on  vortex 
shedding  frequency  was  also  examined. 

An  investigation  on  the  control  of  vortex  shedding  was  also 
done.  A  small-diameter  circular  cylinder  was  maneuvered  parallel 
to  the  main  cylinder  in  the  near  wake  region.  At  low  Reynolds 
numbers  the  control  cylinder  was  able  to  suppress  the  formation 
of  the  vortex  street.  At  a  higher  Reynolds  numbers  the  effect  was 
to  suppress  the  formation  of  three-dimensional  vortex  structures. 
This  control  study  determined  that  actions  in  the  shear  layer  close 
to  the  cylinder  were  responsible  for  the  formation  of  the  vortex 
street.  Furthermore,  the  changes  in  the  wake  shear  layer,  brought 
about  by  the  control  cylinder,  could  be  correlated  to  alterations  in 
the  pressure  field  at  the  cylinder  itself. 
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Overview 


Crystallization  in  polymers  is  a  gradual  process.  In  some  polymers  the  completion  of  this 
process  can  take  several  years.  This  process  can  normally  be  accelerated  by  subjecting  the 
material  to  strain.  In  all  cases  the  characteristics  of  the  resulting  polymer  is  a  function 
of  the  process  under  which  the  crystallization  takes  place. 

To  characterize  this  process  we  have  developed  a  mathematical  model  which  is  pre¬ 
sented  in  chapter..^.  This  model  is  based  on  an  averaging  procedure  that  evaluates  the 
constitutive  equation  for  the  Cauchy  stress  from  averaging  the  effective  response  of  the 
amorphous  (not  crystallized)  part  of  the  material  and  an  array  of  continuously  developing 

crystals. 

In  chapter  3  the  proposed  model  has  been  further  developed  to  represent  those  mate¬ 
rials  which  the  change  in  volume  is  only  related  to  crystallization.  The  model  has  been 

fit  to  results  for  the  crystallization  of  natural  rubber  and  the  predictions  of  the  resulting 
model  for  homogeneous  deformations  has  been  presented.  Included  in  these  results  are 
the  definitions  of  elastic  modulus,  shear  modulus,  and  poisons  ratio  for  a  crystallizing  ma¬ 
terial  which  is  under  large  deformations.  In  addition  to  evaluating  how  these  quantities 
change  in  the  process  of  unconstrained  crystallization  and  crystallization  under  uniaxial 
extension,  the  residual  strain  is  evaluated.  This  part  of  the  work  was  prepared  and  sent 
for  publication. 

The  anomalous  behavior  seen  by  Rosen,  Kolsky,  and  Pipkin  [2]  and  Kolsky  and  Pipkin 
[l]  for  a  torsional  oscillator  with  two  cylindrical  torsional  springs  made  from  natural 
rubber  is  studied  in  chapter  4.  This  study  is  done  in  terms  of  the  proposed  model  and 
it  is  shown  that  the  model  proposed  for  crystallization  can  predict  the  behavior  of  the 
oscillator. 

The  behavior  of  polyethylene  is  studied  in  chapter  5.  An  overview  of  the  work  done 
on  the  mechanical  behavior  of  polyethylene  is  presented.  The  results  from  a  set  of  mod¬ 
ified  stress  relaxation  tests  done  at  10%,  30%,  and  60%  strain  at  room  temperature  are 
presented.  The  modification  to  the  stress  relaxation  was  the  addition  of  small  jumps 
on  the  initial  strain  which  were  used  to  evaluate  the  elastic  modulus  at  fixed  intervals. 
The  results  from  similar  tests  done  at  30%  strain  and  at  35°C,  45°C,  and  55°C  are  also 
presented. 

Concluding  remarks  are  given  in  chapter  6. 
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PROJECT  SUMMARY 


The  metal-fatigue  phenomenon  is  poorly  understood,  partially  because  the 
dislocation  dynamics  of  reversed  deformation  have  not  been  characterized.  This  is 
because  there  has  been,  in  the  past,  only  limited  success  with  direct  observation  of 
dislocations  during  cyclic  deformation,  such  as  with  in-situ  cyclic  or  reversed  plastic 
deformation  tests  in  the  transmission  electron  microscope  (TEM).  This  is  due  to 
problems  associated  with  buckling  of  thin  foils  when  an  applied  tensile  or  shear  stress 
is  reversed.  Analysis  shows  that  dislocation  movement  can  be  reversed  by  tensile 
stressing  in  alternating  perpendicular  directions  (i.e.,  90*  rotations  of  a  tensile  stress  in 
"X"  and  "Y"  directions),  and  buckling  of  the  foil  can  be  avoided.  A  new  specimen 
holder  was  designed  by  the  principal  investigator  for  performing  such  X-Y  mechanical 
tests  in-situ  in  the  TEM.  Recent  tests  supported  by  one-year  funding  by  an  AFOSR- 
Engineering  Foundation  Faculty  Research  Initiation  Grant  and  a  NATO  Collaborative 
Research  grant  have  provided  significant  preliminary  success.  With  some  experimen¬ 
tal  modifications,  it  is  believed  that  a  superior  series  of  in-situ  reversed  deformation 
and  fatigue  experiments  in  the  TEM  can  finally^be  performed.  These  experiments 
would  provide  new  insights  into  reversed  deformation  phenomena,  such  as  the 
Bauschinger  effect  and  fatigue. 
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1.  Introduction  and  background 

Wake  flows  behind  a  variety  of  bodies  ranging  from  circular  cylinders  to  square  rods  and 
screens  to  airfoils  have  been  studied  extensively  (Strouhal  (1S7S),  Rayleigh  (1S79),  Karman 
(1911),  Roshko  (1954),  Sato  and  Kuriki  (1961)  and  others).  Studies  in  the  turbulent  cylinder 
far  wakes  by  Townsend  (1947,  1949)  led  to  the  early  ideas  of  self-preservation  and  Reynolds 
number  independence  (Townsend  1956).  Townsend  postulated  that,  sufficiently  far  down¬ 
stream  of  the  cylinder,  an  asymptotic  self-preserving  state  is  achieved  for  which  the  flow 
can  be  described  by  a  single  velocity  scale  uq  and  a  single  length  scale  L0  (see  figure  1  for 
nomenclature). 

It  has  been  well  established  (Uberoi  Freymuth  1969,  Yamada  et  al.  19S0)  that  far  wake 
flows  are  self-preserving.  It  can  be  shown  that  (see  for  example,  Wygnanski  et  al.  19S6)  the 
nondimensional  velocity  and  length  scales  vary  as 


where,  the  constants  A  and  B  indicate  the  spreading  rate  of  the  wake,  and  x0  is  the  virtual 
origin.  These  constants  would  have  to  be  identical  for  all  wakes  if  the  self- preserving  state 
were  universal. 

Marasli  (1983)  and  Wygnanski  et  al.  (1986)  conducted  experiments  to  determine  whether 
a  universal  self-preserving  wake  exists.  For  this  purpose  a  set  of  experiments  were  designed  to 
investigate  the  wakes  of  porous  screens  of  various  solidities.  A  solid  strip  was  also  considered 
as  a  limiting  case  of  100%  solidity.  The  width  of  the  screens  were  varied  to  keep  the  mo¬ 
mentum  thickness,  6 ,  (hence  the  drag)  constant  at  approximately  2.5  mm.  The  free  stream 
velocity  was  14.5  m/sec.  The  results  are  depicted  in  figure  2  which  shows  the  evolution  of 
the  centerline  velocity  deficit  uo  with  the  nondimensional  streamwise  coordinate.  Solid  lines 
are  linear  least  square  fits  to  the  points  in  the  far  wake  where  self-preservation  is  achieved. 
The  self-preservation  of  each  wake  was  established  from  the  turbulence  intensities,  which  did 
not  vary  with  x  when  normalized  by  the  local  velocity  and  length  scales..  It  is  clear  that 
the  each  of  these  wakes  develop  differently,  indicating  that  the  self-preserving  states  are  not 
universal.  This  is  attributed  to  the  different  initial  conditions  and  the  evolution  of  large  scale 
coherent  structures  in  the  wake.  The  slope  of  each  line  in  figure  2  increases  with  increasing 
screen  solidity.  However,  the  wake  of  the  solid  strip  spreads  the  slowest,  in  contrary  to  the 
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expected  order.  This  behavior  is  quite  puzzling  and  has  not  yet  been  resolved. 

One  possible  explanation  for  this  unexpected  behavior  may  be  given  when  the  near  wakes 
of  these  bodies  are  considered.  Castro  (1971)  reported  that  porous  bodies  with  solidity 
greater  than  80%,  shed  vortices.  If  we  use  his  observation  as  a  gage  and  assume  that  there 
was  no  vortex  shedding  from  the  screens  used  in  our  experiments,  (the  least  porous  screen 
used  had  70%  solidity,)  this  leaves  the  solid  strip  as  the  only  shedding  body.  Monkewitz 
&  Nguyen  (1987)  and  Triantafyllou  et  al.  (1986)  report  predictions  of  the  Karman  vortex 
shedding  frequency  behind  circular  cylinders  via  linear  stability  analysis  of  the  near  wake 
velocity  profile.  Profiles  with  large  velocity  deficits  are  absolutely  unstable  (see  Huerre  and 
Monkewitz  1985)  where  an  impulsively  generated  small  pulse  grows  in  place  and  the  system 
acts  as  a  self-driven  oscillator.  If  the  pulse  is  convected  away,  then  the  profile  is  said  to 
be  convectively  unstable,  which  is  the  case  for  smaller  deficit  wakes.  There  is  encouraging 
evidence  that  shedding  is  a  product  of  absolute  instability  (see  Monkewitz  &  Nguyen  1987). 
It  is  not  clear  why  vortex  shedding  should  affect  the  far  wake  behavior,  since  Karman  type 
vortices  are  known  to  decay  exponentially  with  x.  Nevertheless,  the  possible  contrast  between 
the  far  fields  of  absolutely  unstable  wakes  versus  convectively  unstable  wakes  may  help  clarify 
the  puzzle  stated  earlier. 

An  important  question  remains  unanswered  about  the  relationship  between  the  near  and 
far  wakes.  Can  the  far  wake  of  a  body  be  predicted  from  its  near  field  characteristics?  We 
have  some  evidence  that  (Wygnanski  et  al.  19S6)  the  spreading  rate  of  the  wake  in  the  far 
field  is  related  to  the  near  field  fluctuation  intensity.  Figure  3  shows  the  streamwise  evolution 
of  the  maximum  turbulence  intensities  for  the  screens.  Notice  that  the  70%  screen  has  the 
highest  initial  value,  followed  by  the  45%  and  30%  solidity  screens.  Unfortunately,  the  near 
field  data  for  the  solid  strip  is  not  available,  since  at  the  time  the  main  focus  of  the  research 
was  the  fax  wake.  Nevertheless,  it  seems  like  high  near  field  turbulence  intensity  translates  to 
faster  spreading  of  the  wake,  which  we  proposed  to  establish  for  a  variety  of  wake  generators 
and  Reynolds  numbers. 

It  is  very  promising  that  the  slope  of  the  straight  lines  shown  in  figure  2  can  be  estimated 
from  one  measurement  in  the  near  wake  of  the  generator.  The  prediction  of  the  far  mean 
field  then  would  be  complete  if  the  virtual  origin,  x0,  could  be  obtained  from  the  near  wake 
data.  Recall  that  x0  is  the  x-intcrcept  of  the  straight  line  going  through  the  far  wake  data. 
Wygnanski  'et  al.  (1986)  mention  that  xq  is  used  to  account  for  Reynolds  number  effects. 
However  the  available  data  is  inconclusive  in  any  respect.  The  nondimensional  quantity 
Xo/20  for  the  screen  data  presented  here  are  58,  27,  —21  and  —64  in  increasing  order  of 
solidity,  which  does  not  correlate  with  the  spreading  rates.  We  do  not  know  whether  the 
70%  screen  was  shedding  vortices,  and  earlier  we  discussed  the  possibility  of  its  being  a 
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non-shedding  generator  based  on  the  observations  of  Castro  (1971).  However,  it  it  were 
shedding,  (which  is  quite  possible  since  70%  is  very  close  to  Castro's  threshold  value  of  S0%) 
then  a  negative  xo  would  be  indicating  a  shedding  body.  Marasli  (1983)  reports  negative 
virtual  origins  for  a  circular  cylinder  and  a  roughened  airfoil  (also  known  to  be  shedding) 
and  positive  Xo  for  a  nonshedding  flat  plate.  Then,  the  question  is:  how  does  the  shedding 
characteristics  of  a  wake  relate  to  the  translational  shift  in  its  spreading?  We  claim  that  the 
stability  characteristics  of  the  near  field  fluctuations  are  responsible  from  this  translation 
effect.  Marasli  (1991)  presented  a  discussion  describing  the  basis  for  this  claim. 

Several  important  questions  remain  unanswered  about  the  relation  between  the  near  and 
far  wake  behaviors. 

•  Can  the  spreading  of  the  far  wake  be  predicted  from  the  near  wake  characteristics? 

•  From  far  field  measurements  what  can  we  say  about  the  initial  conditions?  What  can 
we  say  about  the  wake  generator  itself? 

•  What  is  the  physical  meaning  of  the  virtual  origin?  Can  it  be  predicted  from  near 
wake  measurements  only? 

•  What  is  the  effect  of  vortex  shedding  on  the  development  of  the  wake?  Do  absolutely 
unstable  wakes  have  different  far  field  structures  than  convectively  unstable  wakes? 

We  proposed  to  conduct  experiments  in  the  near  and  far  wakes  of  various  wake  generators 
in  order  to  substantiate  the  claims  made  in  the  previous  section  regarding  these  questions. 
We  would  like  to  establish  a  quantitative  link  between  the  initial  conditions  and  the  far  wake 
behavior.  For  this  purpose  we  proposed  a  systematic  investigation  of  wakes  behind  screens, 
flat  plates,  circular  cylinders  etc.  and  correlate  the  wake  spreading  and  the  virtual  origin  to 
the  possible  vortex  shedding  or  the  amplitude  of  the  disturbances  in  the  near  field. 

During  the  initial  stages  of  the  grant  period,  our  efforts  were  concentrated  on  making  the 
wake  measurements  feasible  in  the  University  of  Maryland  Turbulence  Research  Laboratory. 
The  results  of  the  feasibility  studies  were  presented  in  the  progress  report  (Marasli  1991). 
We  also  spent  some  effort  on  the  development  of  a  special  probe  which  would  enable  us  to 
measure  the  vorticity  field  in  the  wake.  Our  initial  attempt  regarding  the  vorticity  probe 
development  was  also  presented  in  Marasli  (1991).  In  this  final  report  we  shall  present  the 
results  of  the  studies  done  on  the  screen  wakes  and  our  more  recent  efforts  on  the  vorticity 
probe  development. 
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2.  Wake  studies  of  screens  with  various  solidities 

Near  wake  studies  have  shown  that  there  are  two  main  types  of  wake  generators:  ones  that 
shed  vortices  into  the  flow  and  ones  that  do  not.  Immediate  flow  fields  behind  these  two  types 
of  wake  generators  are  radically  different  and  lead  to  different  types  of  flow  development  far 
from  the  body.  Unsteady  forces  induced  are  vastly  different  and  at  times  might  lead  to 
resonant  vibrations  of  the  wake  generator. 

We  conducted  experiments  in  the  University  of  Maryland  Turbulence  Research  Labora¬ 
tory  to  study  some  of  the  details  of  the  wakes  behind  screen  with  various  solidities  from  30% 
to  100%.  The  widths  of  the  screens  were  varied  to  keep  the  momentum  thickness  (hence  the 
drag)  constant.  Wakes  were  visualized  using  the  smoke-wire  technique.  Snapshots  of  the 
flow  behind  some  of  the  screens  are  shown  in  figure  4.  The  solidities  of  the  screens  increase 
from  top  to  bottom  and  they  are  30%,  56%,  71%  and  100%  respectively.  The  free  stream 
velocity  was  3  m/s.  It  appears  that  the  wake  width  increases  with  solidity  except  for  the 
100%  case.  This  observation  is  in  accord  with  our  previous  measurements  and  the  reason 
for  this  peculiar  behavior  is  unknown. 

The  near  wakes  of  these  various  generators  are  quite  different.  The  30%  screen  creates 
two  mixing  layers,  one  on  each  side,  which  Initially  evolve  independently  from  one  another 
and  merge  further  downstream  to  form  a  wake.  Shear  layer  instabilities  which  develop  into 
distinct  roll-ups  characterize  this  wake.  The  56%  screen  presents  a  larger  obstacle  to  the 
flow,  hence  the  flow  is  fully  turbulent  immediately  following  the  generator.  Some  traces  of 
organized  motion  is  evident.  The  71%  screen  creates  a  wake  which  is  quite  similar  to  that  of 
a  circular  cylinder.  Vortices  are  shed  from  the  body  at  a  certain  frequency  which  does  not 
change  streamwise  distance.  The  100%  solidity  screen  does  not  show  much  structure  except 
very  close  to  the  body  where  hints  of  von  Karman  shedding  are  present. 

J  ' 

Hot-wire  measurements  of  the  wakes  of  these  screens  are  presented  in  figure  5.  The  quan¬ 
tity  LqUco/uqO,  which  we  shall  call  the  wake  growth  parameter  (Marasli  1989,  Marasli  et  al. 
1991,  1992),  represents  the  spreading  of  the  wake,  and  is  plotted  versus  the  nondimensional 
streamwise  coordinate  x/Q.  The  growth  parameter  can  be  shown  to  vary  iinearly  for  the 
self-preserving  far  wake  ( x/d  >  600.  In  this  phase  of  the  present  project  we  concentrated  on 
the  near  wake,  hence  the  growth  parameter  need  not  vary  linearly  with  x/0.  The  variation 
of  the  growth  parameter  for  the  30%  and  56%  screens  is  shown  in  figure  5a,  while  the  growth 
parameter  for  the  71%  and  100%  screens  is  shown  in  figure  5b.  The  initial  mixing  layer 
region  for  the  30%  screen  is  quite  evident  where  the  growth  parameter  is  constant.  Clearly, 
the  56%  screen  grows  much  faster  thatn  the  30%  screen.  Although  initially  there  is  some 
nonmonotonous  variation  in  the  71%  data,  from  the  last  three  points  it  is  clear  that  the  wake 
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spreading  rate  is  fastest  for  this  case. 

In  order  to  gain  some  insight  to  the  transition  process  from  shedding  to  non-shedding 
wakes,  we  tried  to  determine  the  exact  solidity  where  shedding  starts  to  occur.  In  our  efforts 
we  determined  that  the  transition  occurs  somewhere  between  64%  -  67%  solidity.  Our  next 
goal  was  to  determine  whether  this  transition  is  sudden  or  smooth.  For  this  purpose,  smoke 
photograps  of  the  flow  field  generated  by  the  identical  generator  was  taken  repeatedly.  Figure 
6  shows  photos  of  the  67%  screen-wake  at  a  free  stream  velocity  of  5  m/s  taken  at  various 
instances  of  time.  The  flow  is  dominated  by  von  Karman  type  vortices.  Different  photos 
show  similar  flow  structure  where  the  wavelength  of  these  vortices  are  constant  from  one 
realization  to  the  next.  A  typical  power  spectrum  of  the  streamwise  velocity  fluctuations 
recorded  by  a  hot-wire  in  this  wake  show  a  sharp  peak  at  the  shedding  frequency  (figure  7). 
Short-time  power  spectra  have  also  indicated  that  the  shedding  frequency  does  not  change 
with  time.  Figure  8  shows  photos  of  the  64%  screen-wake.  Clearly  there  is  no  von  Karman 
vortex  shedding  originating  from  the  body.  Organized  structures  seem  to  appear  in  patches 
at  various  streamwise  locations.  The  wavelengths  of  these  structures  are  slightly  different 
in  each  photo.  A  typical  power  spectrum  of  the  streamwise  velocity  fluctuations  recorded 
by  a  hot-wire  in  this  wake  shows  a  broad  band  peak  (figure  9).  Our  initial  attempts  to 
analyse  the  short-time  spectra  has  indicated  that  not  only  these  structures  appear  in  patches, 
their  frequency  vary  with  time  as  hinted  by  the  photos.  This  analysis  however  showed  the 
shortcomings  of  the  Fourier  Transform  as  the  results  were  extremely  difficult  to  interpret. 

The  spectral  peaks  for  the  screens  are  presented  in  figure  10.  The  Strouhal  number, 
Stg  =  fO/Uao  for  the  nonshedding  screens  are  shown  in  figure  10a  versus  x/6.  Notice  that 
the  Stg  decreases  with  downstream  distance,  which  is  the  signature  of  convective  instability. 
Stg  for  the  shedding  screens  is  shown  in  figure  10b.  In  this  case  the  shedding  frequency  and 
the  free  stream  velocities  are  constant  in  each,  case,  thus  the  streamwise  variation  of  Stg  is 
entirely  due  to  uncertainties  in  the  estimation  of  the  momentum  thickness.  Indeed  when  the 
Strouhal  number  based  on  the  screen  width  W,  Stw  =  /VF/C/oo*  is  plotted  we  see  that  it 
remains  constant  in  x  (see  figure  10c.) 

In  conclusion,  several  important  questions  regarding  these  wake  flows  remain  unanswered. 
Recently,  it  has  been  postulated  that  (Huerre  and  Monkewitz  1990)  the  two  distinct  wake 
behaviors  can  be  explained  via  two  different  classes  of  stability,  namely  absolute  and  convec¬ 
tive.  The  vortex  shedding  phenomenon  is  believed  to  be  the  result  of  an  absolute  instability 
originating  near  the  wake  generator  and  causing  the  wake  behave  like  a  self-sustained  oscil¬ 
lator.  By  contrast,  the  non-shedding  wakes  are  convectively  unstable  and  are  characterized 
by  the  absence  of  a  single  sustained  frequency. 


55 


The  relationship  between  the  near  and  far  fields  generated  by  shedding  and  non-shedding 
wake  generators  need  to  be  studied  in  some  detail.  In  order  to  study  the  details  of  the 
transition  from  a  non-shedding  to  shedding  wake,  the  transient  signature  of  the  transition 
has  to  be  recorded.  This  can  be  made  possible  by  looking  at  the  wake  of  two  low  solidity 
screens  placed  one  after  another.  If  the  meshes  are  aligned,  the  composite  solidity  will  be 
low  and  there  won’t  be  any  shedding.  If  the  meshes  are  staggered  the  solidity  will  increase 
and  shedding  will  occur  provided  that  the  composite  solidity  exceeds  the  critical  value  which 
is  estimated  to  be  somewhere  between  64%  and  67%.  The  transient  signal  can  be  recorded 
by  syncronizing  the  screen  movement  with  the  data  acquisition. 

The  short-time  spectral  analysis  of  the  intermittent,  nearly  periodic  signals  via  the  Fourier 
Transform  has  been  observed  to  be  inadequate.  Owing  to  the  uncertainty  of  the  duration 
of  the  intermittencies  the  results  are  often  unreliable.  In  order  to  overcome  this  difficulty 
we  propose  to  apply  the  Wavelet  Transform  technique  (Farge  1992).  The  Wavelet  Trans¬ 
form  is  designed  to  overcome  the  difficulties  encountered  by  the  Fourier  Transform  and  has 
been  successfully  implemented  in  speech  analysis  and  image  processing.  Application  of  the 
Wavelet  Transform  to  turbulent  flows  is  fairly  new  and  is  expected  to  open  new  avenues  in 
the  understanding  of  turbulence. 

3.  The  vorticity  probe 

3.1.  Probe  description  and  calibration  procedure 

Vorticity,  which  is  defined  as  the  curl  of  the  velocity  field  is  known  to  be  the  most 
important  quantity  related  to  the  structure  of  a  turbulent  flow  field.  However,  since  the 
vorticity  involves  the  gradients  of  the  velocity  field,  until  recently  its  measurement  has  been 
limited  to  very  few  simple  cases. 

Vukoslavcevic  et  al.  (1991)  describe  a  miniature  probe  with  nine  hot-wire  sensors  which 
is  capable  of  measuring  the  velocity  and  vorticity  vectors  with  a  spatial  resolution  of  about 
six  Kolmogorov  microscales  in  a  a  turbulent  boundary  layer. 

The  orientation  of  the  sensors  and  the  calibration  procedure  was  believed  to  hinder  the 
performance  of  the  Vukoslavcevic  probe.  Hence,  a  new  circular  design  and  calibration  tech¬ 
nique  was  developed  by  Marasli  (1991).  This  new  probe  was  designed  to  measure  nine  un¬ 
knowns  which  are  the  three  components  of  the  instantaneous  velocity  vector  u  =  (u,  v ,  in)  and 
the  six  instantaneous  velocity  gradients  du/dy,  du/dz,  dv/dy,  dv/dz,  dw/dy  and  dw/dz. 
The  gradients  in  the  streamwise  direction  are  then  obtained  via  Taylor’s  hypothesis  (Piomelli 
et  al.  1989),  thus  enabling  the  estimation  of  all  three  components  of  the  instantaneous  vor¬ 
ticity  vector.  The  response  equations  obtained  from  the  nine  sensors  would  provide  the 
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information  required  the  solve  for  the  nine  unknowns.  The  testing  of  this  newly  circular 
probe  showed  the  basic  weakness  of  any  nine-sensor  probe:  since  the  hot-wire  response 
equations  are  nonlinear,  the  solution  of  a  9x9  system  is  not  unique.  Which  means  that  for 
a  given  set  of  anemometer  voltages,  there  correspond  more  than  one  solution,  which  is  not 
acceptible  physically.  Extra  information  is  needed  to  resolve  the  nonuniqueness. 

A  yet  newer  probe  which  consists  of  twelve  sensors  was  developed.  The  twelve  sensors 
provide  twelve  equations  to  solve  for  nine  unknowns.  The  redundant  information  is  used  to 
resolve  the  nonuniqueness  problem.  The  sketch  of  the  new  probe  is  shown  in  figure  11.  The 
twelve  sensors  are  arranged  in  three  arrays  each  consisting  of  four  wires.  The  centers  of  the 
arrays  are  located  at  the  corners  of  an  equilateral  triangle  which  measures  1.32  mm  on  each 
side.  The  prongs  on  which  the  wires  are  welded  are  separated  by  0.44  mm.  The  24  tungsten 
prongs  are  tapered  down  to  30  fi  diameter  at  the  tip.  The  sensors  are  also  made  of  tungsten 
and  are  2.5  fi  in  diameter. 

The  quality  of  any  hot-wire  measurement  technique  starts  with  the  quality  of  the  calibra¬ 
tion.  Regardless  of  the  technique  used,  the  calibration  of  hot-wire  probes  deteriorate  with 
time  due  to  various  reasons  such  as  temperature  variations,  drift  in  electronics,  impurities  in 
the  flow,  etc.  Therefore,  speedy  calibration  and  measurement  are  essential.  The  technique 
described  here  enables  the  calibration  of  a  multi-sensor  probe  (non-zero  pitch-yaw  combi¬ 
nations  included)  in  less  than  5  minutes  with  a  computer  controlled  traversing  mechanism. 
The  short  duration  of  calibration  period  directly  translates  to  improved  measurement  qual¬ 
ity.  Following  Marasli  and  Nguyen  (1991),  we  start  with  the  effective  cooling  velocity  which 
is  given  by  Jorgensen’s  (1971)  response  equation: 

u;„  =  u%+ tvT  +  hVB  (3) 


where  ujy,  Uj  and  u2g  are  the  normal,  tangential  and  binormal  components  of  the  velocity 
with  respect  to  the  sensor  (see  figure  lib)  and  k  and  h  are  the  tangential  and  binormal 
cooling  coefficients  respectively.  A  fourth  order  polynomial  is  used  to  relate  the  effective 
cooling  velocity  to  the  anemometer  bridge  voltage  in  the  following  form: 

U2ff  =  P{E)  =  A0  +  AXE  +  A2E2  +  A3E3  +  AaE\  (4) 

For  a  slanted,  inclined  wire 


nxu  +  n2v  -f  n3w, 

(5) 

t\u  +  t2v  +  t3w , 

(6) 
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uB  =  b\u  +  b2v  +  b3w.  (7) 

where  n,-,  £,•  and  6,  (i  =  1-3)  are  known  coefficients  that  depend  on  the  orientation  of  the 
wire.  Substitution  of  (5-7)  into  (4)  yields: 

Uljj  =  A\U 2  +  A2V2  +  A3w2  +  A4uv  +  A$uw  +  Aevw  =  P(E )  (8) 

For  each  of  the  twelve  sensors  in  the  probe  we  end  up  with  an  equation  such  as  (8).  Since 
we  have  nine  unknowns,  (three  velocity  components  and  six  gradients,)  the  twelve  equations 
provide  extra  information  which  can  be  utilized  to  resolve  the  nonuniqueness  problem.  We 
choose  to  solve  the  equations  in  a  least-square  sense,  whicli  provides  unique  solutions  within 
the  calibration  domain.  For  convenience  we  write  (S)  as 

u2  =  Pj{Ej)  +  K\v2  +  K2w2  +  K3u.jVj  +  K^ujWj  +  K$VjWj  (9) 

for  j  =  1-12  sensors.  Notice  that  we  have  lumped  the  known  coefficients  n,,  t,  and  6,  into 
the  unknown  calibration  coefficents  A'tj  to  be  determined  during  calibration.  Thus  errors 
due  to  uncertainties  in  the  slant  and  roll  angles  are  eliminated,  which  is  one  major  advantage 
of  this  calibration  procedure.  When  a  multi-sensor  probe  is  used  in  a  non-uniform  flow,  the 
velocity  vector  seen  by  each  sensor  will  be  different.  Following  Vukoslavcevic  et  al.  (1991) 
we  expand  the  velocity  measured  by  each  sensor  in  a  Taylor  series  around  the  centroid  of 
the  probe.  Hence  we  let 


Vj  =  V  +  c 


du  du 

%  +  diTz+" 

(10) 

dv  dv 

’aZ+d’Tz+"' 

(ii) 

dw  ,  dw 
+  diTz  + ' " 

(12) 

where  the  right-hand-side  is  evaluated  at  the  probe  centroid.  The  constants  c,  and  dj  are 

the  vertical  and  horizontal  distances  from  the  centroid  of  each  wire  to  the  centroid  of  the 
probe  respectively.  Substitution  of  (10-12)  into  (9)  yields  12  nonlinear  algebraic  equations 
with  the  9  unknowns-  - 


du  du  dv  dv  dtv  dw 

U,  V ,  W ,  “  y  1  r\  1  "a  and  “o 

ay  uz  ay  az  ay  oz 


(13) 


evaluated  at  the  centroid  of  the  probe,  which  can  be  solved  in  a  least-square  sense  using 
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standard  numerical  techniques  such  as  Newton's  method.  The  nonlinear  equations  are: 


du  du  dv  dv  dwdw 

/,(u, u, 07’  ~fz'  j;- P>'c>'  ~  0 


where 


Pj  =  Aoj  +  A\jEj  +  AijE2j  -f  AijEj  +  A^jE* 


and 


,  n  2  dU 

fj  =  -Pj  +  U  +  2  CjU—  -f  2dj  u  — 

dy  A~ 


-A'i  i 
-A% 
-A3j 

-/Csi 


<9u 


&r 

<9u 


„  +2cjt,_+2<i)„-J 

iu2  +  2Cjty—  +  2djW~^ 

ay  oz 

(  dv  du)  (  dv  dus 

a.  A  ( „ 

dz 


(  dw  du)  (  dw  du'' 

uw  +  c’  (“  ej  +  w  Ty  )  + <H“' Tz  +  wTz, 

dv 


(  dw  dv)  (  dw 

TO+Cjit)_,.„_i+<i,.it,_+K,& 


) 


(14) 


(15) 


(16) 


(17) 


The  unknown  constants  Aij  and  Kij  are  obtained  during  calibration  by  subjecting  the 
probe  to  uniform  flow  (Q)  coming  at  various  pitch  (a)  and  yaw  ((3)  angles.  The  calibration 
velocity  components  (u,v,w)  are  then  given  by  u  =  Q  cos  a  cos/?,  v  =  Q  sin  a  cos/?  and  w  ~ 
Q sin/?.  For  calibration  the  probe  is  mounted  on  the  traversing  mechanism  of  a  calibration 

W  * 

jet  which  keeps  the  sensing  volume  of  the  probe  at  the  same  location  as  the  pitch  and  yaw 
angles  are  varied.  Typically,  a  combination  of  3  yaw  and  3  pitch  angles  are  used  between 
±20°  at  two  different  speeds,  thus  giving  a  total  of  18  calibration  points.  A  typical  calibration 
curve  for  one  sensor  is  shown  in  figure  12  where  the  calibration  data  points  (indicated  by 
symbols)  are  plotted  versus  the  bridge  voltage.  For  the  symbols,  the  ordinate  is  U2jj,  where 
the  velocity  vector  u  is  known  from  the  pitch  and  yaw  angles,  and  the  constants  Ki  have 
been  obtained  during  calibration.  The  solid  curve  is  a  fourth  order  polynomial  fit  which 
describes  the  relationship  between  the  effective  cooling  velocity  and  the  bridge  voltage.  The 
fact  that  all  points  collapse  on  one  curve  is  an  indication  that  the  assumed  cooling  law  is 
satisfactory.  Figure  12  displays  the  calibration  curve  for  one  sensor  only.  All  twelve  sensors 
of  the  probe  exhibit  similar  calibration  properties.  The  quality  of  the  calibration  can  be 
checked  by  trying  to  recover  the  velocities  used  in  the  process.  In  addition,  since  the  probe 


59 


is  calibrated  in  a  uniform  flow,  there  should  be  no  gradients  in  the  calibration  data.  A  good 
calibration  is  one  which  recovers  the  velocities  within  1%  and  gives  spurious  vorticities  that 
are  less  than  20  sec-1.  The  sample  data  shown  here  fall  in  that  category. 

Note  that  the  present  calibration  scheme  can  easily  be  applied  to  probes  with  different 
geometries.  The  only  geometrical  information  needed  is  the  distance  from  the  center  of  each 
wire  to  the  centroid  of  the  frontal  sensing  area  of  the  probe.  The  number  of  sensors  need 
not  be  fixed  at  twelve. 

3.2.  Vorticity  measurements  in  the  wake  of  a  circular  cylinder 

The  performance  of  the  twelve-sensor  vorticity  probe  was  tested  in  the  wake  of  a  6.35  mm 
diameter  circular  cylinder.  The  probe  was  place  downstream  of  the  cylinder  at  the  location 
x/d  —  30.  The  free  stream  velocity  was  =  5  m/s,  resulting  in  a  Reynolds  number  of 
2000  based  on  the  diameter. 

The  mean  velocity  profile  is  shown  in  figure  13.  The  ordinate  is  the  vertical  coordinate 
normalized  by  the  half-width  of  the  wake,  Lq ,  which  was  S.5  mm.  The  abscissa  is  the  velocity 
deficit  normalized  by  the  centerline  value  u0  =  0.S2  m/s.  The  momentum  thickness,  6,  which 
is  proportional  to  the  drag  of  the  cylinder  was  2.S  mm.  In  figure  13  the  symbols  represent  the 
measurements  and  the  solid  curve  corresponds  to  an  analytical  curve  fit  tho  the  experimental 
data.  The  profile  is  quite  symmetric  around  the  centerline  as  it  should  be. 

The  turbulence  intensities  of  the  u  (streamwise  or  x-direction),  v  (cross-stream  or  in¬ 
direction)  and  w  (spanwise  or  z-direction)  velocity  components  are  shown  in  figure  14.  The 
v  component  has  the  most  intense  fluctuations,  which  is  due  to  the  von  Karman  vortex 
shedding  present  in  the  flow.  The  Reynolds  shear  stresses  are  presented  in  figure  15.  Notice 

that  the  u'u /  and  v'w'  components  are  negligible  compared  to  the  dominant  uV  component. 
___ _____  *  ' 

Ideally  u'w'  and  v'w'  should  by  identically  zero. 

Figure  16  compares  the  distribution  of  the  mean  velocity  gradient  dU  jdy  measured  di¬ 
rectly  with  the  twelve-sensor  probe  with  the  derivative  of  the  curve  fit  to  the  U  distribution 
shown  in  figure  13.  Notice  the  excellent  agreement  between  the  two,  which  shows  the  capa¬ 
bility  of  this  probe  to  measure  instantaneous  velocity  gradients  accurately,  which  is  essential 
for  vorticity  measurements. 

The  mean  vorticity  components  are  shown  in  figure  17.  Notice  that  fL  component  is  the 
dominant  one  which  results  from  the  mean  velocity  gradient  ion  the  cross-stream  direction. 
The  Rr  and  fly  components  should  ideally  be  zero.  The  extremely  small  nonzero  components 
recorded  is  an  overall  indication  of  the  lack  of  spanwise  uniformity  in  the  flow  and  the  ability 
of  the  probe  to  resolve  the  gradients.  The  variance  of  the  vorticity  fluctuations  are  shown  in 
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figure  IS.  Unlike  the  mean  vorticitv.  all  three  components  of  the  vorticity  fluctuations  are 
of  the  same  order  of  magnitude,  which  is  a  property  of  turbulence  in  general. 

The  power  spectra  of  the  velocity  fluctuations  at  three  different  cross-stream  locations 
are  shown  in  figure  19.  The  sharp  peak  in  the  spectra  corresponds  to  the  von  Karman  vortex 
shedding,  which  occurs  at  160  Hz  or  Strouhal  number,  St  =  fd/Uoo  —  0.2.  Notice  that  the 
shedding  is  most  evident  in  the  u-spectra  and  nonexistent  in  the  tu-spectra,  which  is  an 
indication  of  the  nominal  two-dimensionality  of  the  flow.  The  vorticity  spectra  are  shown 
in  figure  20.  Again  the  shedding  frequency  is  visible  in  the  spanwise  vorticity  component, 
which  corresponds  to  the  two-dimensional  nature  of  the  shedding. 

In  conclusion,  we  consider  these  measurements  in  the  circular  cylinder  wake  as  a  landmark 
success.  To  our  knowledge,  this  is  the  first  documentation  of  the  simultaneous  measurement 
of  all  three  components  of  the  velocity  and  all  three  components  of  the  vorticity  vectors. 
Efforts  toward  the  usage  of  this  probe  in  flows  where  the  three  dimensional  structure  is  more 
eminent  is  underway. 

4.  Publications  resulting  from  present  grant 

Funds  provided  by  the  present  grant  were  crucial  in  the  continuation  of  the  research  projects 
which  resulted  in  the  publication  Marasli  &  Nguyen  (1991).  We  also  submitted  an  abstract 
to  2nd  International  Symposium  on  Engineering  Turbulence  Modelling  and  Measurements 
to  be  held  on  May  31  -  June  2,  1993  in  Florence,  Italy.  Mr.  Phuc  Nguyen  is  doing  his 
Ph.D.  dissertation  on  the  vortical  structure  of  the  wake  and  his  dissertation  work,  upon 
completion,  will  be  submitted  to  an  archival  journal  for  publication.  Mr.  John  Wright,  who 
is  doing  his  M.S.  thesis  on  the  screen  wakes,  will  finish  soon  and  his  work  will  be  submitted 
to  an  archival  journal  for  publication.  The  sponsorship  of  the  Engineering  Foundation  will 
be  gratefully  acknowledged  in  those  publications. 

5.  References 

Castro,  I.  1971  J.  Fluid  Mech.  46,  599. 

FaRGE,  M.  1992  “Wavelet  transforms  and  their  applications  to  turbulence,”  Annu.  Rev. 
Fluid  Mech.,  24,  395. 

HUERRE,  P.  &  MoNKEWITZ,  P.  A.  1985  J.  Fluid  Mech.  159,  151. 

HUERRE,  P.  &:  MoNKEWITZ,  P.  1992  “Local  and  global  instabilities  in  spatially  developing 
flows,”  Annu.  Rev.  Fluid  Mech.,  24,  395. 

Jorgensen,  F.  E.  1971  DISA  Inform.  11,  31. 

KaRMAN,  T.  Von  1911  Gottinger  Nachr.,  math-phys.  Kl.,  509. 


61 


MaRASLI,  B.  1983  Master's  thesis.  Dept,  of  AME,  University  of  Arizona. 

MaRASLI,  B.  1991  The  structure  of  wakes  behind  various  bluff  and  non-bluff  bodies. 

Progress  Report  submitted  to  Engineering  Foundation. 

MaRASLI,  B.  1989  Ph.D.  dissertation,  University  of  Arizona. 

Marasli,  B.,  Champagne,  F.  Sc  Wygnanski,  I.  1991  Phys.  Fluids  A  3,  665. 
Marasli,  B.,  Champagne,  F.  &  Wygnanski,  I.  1992  J.  Fluid  Mech.  235,  511. 
MARASLI,  B.  Sc  Nguyen,  P.  1991  “A  calibration  technique  for  multiple-sensor  hot-wire 
probes,”  Bull.  Am.  Phys.  Soc.  36,  2623. 

Monkewitz,  P.  A.  Sc  Nguyen,  L.  N.  1987  J.  Fluids  Struct,  l,  165. 

PlOMELLl,  U.,  BALINT,  J.-L.  Sc  Wallace,  J.  M.  19S9  Phys.  Fluids  32,  609. 
Rayleigh,  Lord  1879  Phil.  Mag.  ser.  5,  7(42),  149. 

Roshko,  A.  1954  NACA  Rep.  1191. 

Sato,  H.  Sc  Kuriki,  K.  1961  J.  Fluid  Mech.  11,  321. 

STROUHAL,  V.  1S78  Ann.  Phys.  Chem.  5,  217. 

Townsend,  A.  A.  1947  Proc.  R.  Soc.  Lond.  A190,  551. 

Townsend,  A.  A.  1949  Austral.  J.  Sci.  2,  451. 

TOWNSEND,  A.  A.  1956  The  Structure  of  Turbulent  Shear  Flow.  Cambridge  University 
Press. 

Triantafyllou,  G.  S.,  Triantafyllou,  M.  S.  Sc  Chryssostomidis  C.  1986  J.  Fluid 
Mech.  170,  461. 

Uberoi,  M.  S.  Sc  FreYMUTH,  P.  1969  Phys.  Fluids  12,  135T 

Vukoslavcevic,  P.,  Wallace,  J.  M.  Sc  Balint,  J.-L.  1991  J.  Fluid  Mech.  228,  25. 
Wallace,  J.  M.  1983  Bull.  Am.  Phys.  Soc.  28,  1401. 

Wygnanski,  I.  J.,  Champagne,  F.  H.  Sc  Marasli,  B.  1986  J.  Fluid  Mech.  168,  31. 
Yamada,  H.,  Kawata,  Y.,  Osaka,  H.  Sc  Kageyama,  Y.  1980  Tech.  Rep.  Yamaguchi 
University,  Japan,  2(4). 


62 


A  30%  Solidity  screen 
A  45  %  " 

G  70  %  " 


65 


Figure  3.  The  maximum  turbulence  intensities  versus  the  strenmwise  coordinate  for  the  screens. 


FIGURE  5.  Variation  of  the  mean  flow  growth  parameter  with  streamwise  distance  for  the 
screen  wakes. 


Figure  ^  .  Smoke  photographs  behind  the  wake  of  64o  solidity  screen. 

Each  photo  corresponds  to  a  different  instant  in  time. 


70 


Figure 


Ste 


Stg 


e  zee  4ee  see  see 


Stw 


0.Z5 

0.2 

0.1S 

0.1 

0.05 

0 

0  200  400  600  800 


i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 


^  ooo 


_i _ I _ I _ I _ L. 


O  71* 
O  100S 

II  II  I  I  I  I 


x/0 


.  Variation  of  the  Strouhai  number  with  streamwise  distance  for  the  screen  wakes. 
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Typical  calibration  data  for  one  sensor  of  the  vorticity  probe.  The  symbols 
represent  the  measured  data  while  the  solid  curve  corresponds  to  the  polynomial 
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FIGURE  13.  Mean  velocity  profile  measured  by  the  vorticity  probe  in  the  circular  cylinder 
wake. 
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FIGURE  14.  Turbulence  intensities  measured  by  the  vorticity  probe  in  the  circular  cylinder 
wake. 
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FIGURE  15.  Reynolds  shear  stress  profiles  measured  by  the  vorticity  probe  in  the  circular 
cylinder  wake. 
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FIGURE  16.  Mean  flow  gradient  measured  by  the  vorticity  probe  in  the  circular  cylinder 
wake. 
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17.  Mean  vorticity  profiles  measured  by  the  vorticity  probe  in  the  circular  cylinder 
wake. 
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FIGURE  18.  Variance  of  vorticity  fluctuations  measured  by  the  vorticity  probe  in  the  circular 
cylinder  wake. 
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REPORT  SUMMARY 


This  report  documents  the  results  of  a  one  year  project  funded  by  Air  Force  Engineering 
Foundation  Research  Initiation  Grant  RI-B-90-15.  A  general  framework  for  small  strain, 
isothermal,  cyclic  viscoplasticity  of  metals  and  alloys  is  presented.  The  framework  employs 
multiple  nonlinear  kinematic  hardening  rules.  An  approach  for  partitioning  isotropic  harden¬ 
ing  effects  between  the  so-called  “yield  stress”  and  “drag  stress”  in  the  viscoplastic  flow  rule 
and  the  saturation  level  of  the  kinematic  hardening  subvariables  is  presented.  A  bounding 
surface  interpretation,  which  is  useful  for  interpretation  of  experimental  results  is  also  given. 
Materials  and  biaxial,  tension-torsion  experiments  on  thin  walled  tubular  specimens  are  de¬ 
scribed  and  then  used  to  illustrate  a  technique  for  estimating  overstress  from  uninterrupted 
cyclic  tests.  The  results  of  these  estimations  are  interpreted  in  the  context  of  the  general 
framework.  The  technique  is  applied  to  data  obtained  for  Type  304  stainless  steel  tested 
at  room  temperature  and  for  Waspaloy  (fine  and  coarse  grain  microstructures)  tested  at 
650°C.  Determination  of  the  appropriate  number  of  nonlinear  kinematic  hardening  rules  is 
accomplished  by  differentiation  of  backstress  histories  for  proportional  cycles  computed  with 
the  knowledge  of  the  effective  overstress  at  approximately  the  same  effective  inelastic  strain 
rate.  Special  cases  of  the  general  nonlinear  kinematic  hardening  rule  are  considered.  The 
appropriate  distribution  of  isotropic  hardening  among  the  saturation  levels  of  the  nonlinear 
kinematic  hardening  subvaxiables  is  considered  in  light  of  the  presented  results.  Kinematic 
hardening  rates  are  compared  with  experimental  results  and  refinement  is  made  based  on 
this  comparison.  Isotropic  hardening  rules  and  rate  sensitivity  functions  are  formulated  and 
the  complete  model  is  integrated  for  representative  biaxial,  strain-controlled  test  histories 
and  compared  with  experimental  results  for  both  Waspaloy  microstructures. 
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V. 


CONCLUSIONS 


It  can  be  concluded  from  the  results  of  this  study  that  the  implementation 
of  the  unified  state  variable  constitutive  theory  into  the  ADINA  finite  element 
program  yields  important  information  pertaining  to  the  mechanical  behavior  of  the 
fiber/matrix  system.  From  the  preliminary  tensile  tests  of  the  constitutive 
model  in  ADINA,  it  was  found  that  the  model  had  the  capability  of  capturing  such 
material  response  as  the  stress  drop  phenomenon  which  would  not  be  observed  in 
a  bilinear  elastic-plastic  model.  A  reasonable  finite  element  mesh  was  selected 
after  extensive  comparisons  of  a  variety  of  models. 

As  a  result  of  the  research,  insight  has  been  gained  as  to  the  magnitude 

of  the  residual  stress  field  located  at  the  fiber/matrix  interface.  Various 

cooling  profiles  were  investigated  and  shown  to  have  a  limited  effect  on  the 

reduction  of  the  residual  stress  field.  Stress  reductions  achieved  through 

temperature  and  pressure  holds  of  various  times  were  also  seen  to  change  the 

residual  stress  field  very  little.  Stress  ratchetting  was  observed  in  thermally 

cycled  composites.  By  applying  a  constant  radial  pressure  to  the  composite  and 

thermally  cycling  the  system  between  300°C  and  600°C,  a  reduction  in  the  residual 
stress  field  magnitude  was  achieved.  This  thermal-  and  pressure -loading 

combination  shows  the  most  promise  for  the  development  of  composites  with 

residual  stress  fields  below  the  magnitudes  at  which  fiber/matrix  Interface 

damage  is  initiated. 
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ABSTRACT 


An  investigation  was  conducted  to  study  the  stress  corrosion 
cracking  susceptibility  of  superplastically  deformed  Weldalite  049, 
which  was  available  in  a  thermomechanically  processed  form, 
suitable  for  superplastic  forming.  Tensile  test  specimens  machined 
from  sheets  of  0.090  inch  thickness  were  initially  superplastically 
deformed  to  different  strain  levels  at  a  constant  true  strain  rate 
of  2xlO'Vs  and  a  temperature  of  about  500  C  using  an  INSTRON 
superplastic  testing  facility,  available  in  the  Department  of 
Mechanical  Engineering  of  Texas  A&M  University.  Most  specimens  were 
deformed  under  a  hydrostatic  (back)  pressure  of  400  psi,  but  some 
selected  specimens  were  deformed  without  back  pressure.  After 
superplastic  deformation,  most  specimens  were  aged  by  solutionizing 
at  504  C  for  1  hr,  followed  by  precipitation  hardening  for  24  hr  at 
180  C,  but  again  some  selected  specimens  were  not  aged. 

Stress  corrosion  testing  was  then  conducted  on  the 
superplastically  deformed  specimens  at  an  (initial)  slow  strain 
rate  of  2xl06/s,  using  a  CORTEST  Slow  Strain  Rate  Testing  Machine. 
A  test  cell  specially  designed  and  fabricated  by  the  investigator 
was  used  for  this  purpose.  The  tests  were  conducted  using  3.5%  NaCl 
solution  in  deionized  water  as  the  electrolyte.  The  test  potential 
was  mostly  -800  mV,  but  in  selected  cases  the  tests  were  conducted 
at  other  potentials.  For  comparative  purposes,  tests  were  also  made 
on  as-received  samples,  both  in  the  electrolyte  and  in  air. 
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The  results  indicate  that  for  the  aged  samples 

superplastically  tested  under  back  pressure,  the  elongation  at 
-800  mV  is  more  sensitive  to  superplastic  strain  than  the  tensile 
strength,  though  both  the  values  show  a  decreasing  trend  with 
increase  in  the  superplastic  strain.  The  application  of  back 
pressure  during  superplastic  testing  helps  to  improve  the  stress 
corrosion  resistance,  presumably  by  resisting  the  growth  of 
cavities  and  the  activation  of  potential  cavity  nuclei.  Aging 
impi oves  the  inherent  strength  of  the  material  and  is  therefore 
desirable  for  improving  the  resistance  to  stress  corrosion 
susceptibility  in  relation  to  not  aged  samples.  The  aged  samples 
show  increased  susceptibility  for  stress  corrosion  cracking  as  the 
applied  potential  is  moved  towards  the  pitting  potential.  Under 
pitting  conditions,  these  specimens  are  completely  corroded. 
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ABSTRACT 


A  technique  is  presented,  which  renders  computational  the  bounds  on 
the  effective  properties  of  two-dimensional  composites  developed  by 
previous  researchers.  Statistical  integrals  appearing  in  these  bounds 
contain  three-point  correlation  functions  on  material  properties  in  the 
plane  transverse  to  the  fiber  axis.  The  correlation  functions  can  be 
systematically  reexpressed  in  terms  of  the  substructure  geometry,  for 
fiber  reinforced  composites:  matrix  properties,  fiber  size  and  relative 
positioning,  and  fiber  properties.  The  method  presented  here  is 
innovative  in  that  it  allows  for  non-uniform  material  properties  within 
the  fibers  themselves.  The  motivation  for  this  is  to  predict  properties 
of  composites  when  there  is  an  arbitrary  transition  into  the  surrounding 
matrix,  such  as  a  coating,  a  surface  roughness  which  creates  a  smooth 
transition  region,  or  an  incipient  separation  of  the  fiber-matrix 
interface.  The  usual  model  case  of  uniform  fibers  is  treated  first,  and 
then  the  more  general  case  of  non-uniform  fibers. 
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5.0  CONCLUSIONS 


During  this  project,  a  number  of  conformal  antenna  geometries  have  been  studied  as 
a  means  of  overcoming  some  of  the  disadvantages  of  conventional  microstrip  patch 
antennas  such  as  narrow  bandwidth,  restricted  scan  volume,  and  large  echo  area  at  discrete 
frequencies  without  comprising  attractive  features  such  as  low  cost  and  light  weight 
Using  the  Green's  function/Galerkin  technique,  numerical  models  for  single  radiating 
elements  comprising  circular  patches  backed  by  rectangular  cavities  and  circular  patches 
surrounded  by  conducting  super-planes  and  for  infinite  arrays  of  circular  patches  with  each 
patch  backed  by  a  rectangular  cavity  have  been  developed.  Results  from  these  numerical 
models  have  been  compared  to  well-validated  results  for  equivalent  conventional  circular 
patch  cases.  The  new  geometries  exhibit  impedance  characteristics  and  radiation  patterns 
within  the  operating  band  that  do  not  differ  greatly  from  those  of  conventional  circular 
patches.  However,  scattering  characteristics  for  the  new  geometries  differ  significantly 
outside  of  the  operating  band  from  those  of  conventional  circular  patches.  They  possess  an 
extra  degree  of  freedom  which  allows  minimization  of  the  RCS  at  discrete  frequencies 
outside  of  the  operating  band,  which  may  be  useful  in  certain  applications.  Perhaps  the 
most  important  result  of  this  research  has  been  the  confirmation  of  the  fact  that  the  scan 
performance  of  conformal  patch  arrays  can  be  dramatically  improved  by  breaking  the 
substrate  into  discrete  cavities.  Indeed,  cavity-backed  patch  arrays  do  not  suffer  from  the 
scan-bandwidth  tradeoff  inherent  to  conventional  patch  arrays.  In  summary,  it  seems  that 
Geometry  2  is  useful  for  situations  where  it  is  desirable  to  minimize  scattering  at  certain 
frequencies  outside  the  operating  band  of  the  antenna  It  is  better  for  this  situation  than 
Geometry  1  because  the  additional  super-plane  conductor  should  add  very  little  weight  to 
the  antenna  where  as  a  cavity-backing  would  probably  add  substantial  weight.  For  arrays, 
the  dramatic  improvement  in  scan  performance  probably  overcomes  the  drawback  of  the 
additional  weight  of  Geometry  1.  It  is  interesting  to  note  that  other  workers  are  also 
beginning  to  consider  cavity-backed  patch  geometries  [2].  In  [2J,  another  advantage  of  the 
cavity-bae'  1  geometry  is  mentioned:  improved  heat  dissipation  which  is  important  when 
active  devices  are  used.  The  cavity-backing  described  in  [2]  is  circular,  and  in  the  future, 
we  intend  to  extend  our  solution  to  circular  cavity  backings  because  these  may  be  easier  to 
machine. 

Evaluation  of  the  RCS  of  these  various  geometries  requires  substantial  CPU  time. 
We  are  currently  looking  into  various  interpolation  techniques  which  would  allow  the  RCS 
to  be  accurately  interpolated  at  very  fine  frequency  increments  from  the  exact  calculation 
which  is  performed  only  over  a  very  few  frequencies.  One  promising  technique  seems  to 
be  model-based  parameter  estimation  [3]. 

We  are  also  attempting  to  accelerate  the  convergence  of  the  series  expressions  for 
the  matrix  element  zn+j  for  Geometry  1  using  similar  techniques  as  are  used  for 
Geometry  2  in  Section  2.3.  However,  results  to  date  have  not  been  promising. 

Experimental  work  on  tunable  conformal  antennas  to  date  shows  some  promise, 
and  work  is  continuing  in  this  area.  Perhaps  the  most  desirable  achievement  would  be  the 
development  of  an  analytical  model  for  a  slot  backed  by  a  cavity  that  is  partially  filled  with 
ferrite. 
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UMMARY  AND  DISCUSSION 
emulations  have  been  presented  of  TOF  ex- 
nents  across  abrupt  and  graded  amorphous 
[conductor  heterojunctions.  The  features  of 
transients  can  be  associated  with  processes 
h  take  place  in  the  device  during  simulated 
iriments.  If  the  transients  were  more 
plicated,  their  features  could  be  analyzed  by 
ying  internal  variables  of  the  simulator. 

'he  simulations  suggest  that  TOF  expe¬ 
rts  can  make  an  important  contribution  to 
understanding  of  amorphous  semiconductor 
xojunctions.  Experiments  with  abrupt 
:tions  and  with  various  graded  junctions 
d  show  if  an  abrupt  heterojunction  is  really 
ivalent  to  a  graded  junction  of  zero  or  near 
width.  If  not,  the  difference  would  be  due  to 
effect  of  surface  states  at  the  heterojunction, 
experiments  may  also  provide  a  direct  mea- 
t  of  the  band  energy  change  at  the  junction. 

5y  studying  capture  and  emission  rates  as 
i  in  TOF  transients  with  different  device 
knesses,  it  should  also  be  possible  to  find 
e  reliable  values  of  the  capture  coefficients 
the  states  at  different  energies.  In  addition, 
sriments  such  as  these  may  give  information 
tow  shallow  tail  states  contribute  to  transport 
region  where  the  mobility  edge  is  shifting  to 
igher  energy.  This  may  reveal  important 
ts  to  the  applicability  of  the  band  conduction 
[el  to  room  temperature  conduction  in 
rphous  materials. 

)nce  the  experiments  are  performed,  it  will 
>ossible  to  simulate  them  using  the  same 
nimental  parameters.  It  will  then  be 
ible  to  make  direct  comparisons  between  the 
.  and  the  simulation  parameters. 
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Project  Description 

This  project  addresses  the  need  for  a  well-developed  theory  in  the  area  of 
continuous-time,  or  analog,  adaptive  signal  processing.  To  date,  the  adaptive  processing  of 
signals  has  been  conducted  almost  entirely  with  discrete-time,  or  digital,  implementations. 
Although  the  digital  implementation  is  effective  for  a  broad  class  of  signals,  it  has  two 
drawbacks.  Firstly,  many  signals  exist  which  have  bandwidths  beyond  the  capabilities  of 
digital  hardware.  Although  the  speed  of  digital  circuitry  is  increasing  rapidly,  the 
bandwidths  of  signals  in  telecommunications,  radar,  etc.,  can  always  be  expected  to  expand 
just  as  rapidly,  thus  remaining  out  of  the  reach  of  digital  adaptive  filters.  For  the  same 
reasons,  analog  circuitry  will  always  have  bandwidth  capabilities  beyond  those  of  digital 
electronics,  making  analog  implementations  attractive  for  wideband  adaptive  filters. 
Secondly,  for  a  given  application,  analog  circuitry  is  much  more  efficient  than  digital 
circuitry  in  terms  of  complexity,  size,  weight,  and  power.  Thus,  if  performance  is  equivalent 
for  the  two  implementations,  the  analog  implementation  is  usually  preferable.  Demonstra¬ 
tion  of  this  equivalence  of  performance  is  one  of  the  goals  of  this  project,  and  success  in  this 
respect  has  been  achieved. 

Specific  objectives  for  the  project  are  as  follows: 

1)  Comparison  of  discrete-time  (DT)  and  continuous-time  (CT)  adaptive  filters.  The 
first  step  in  this  task  is  the  development  of  algorithms  for  the  adaptation  of  the  CT  filter. 
Following  this  development  is  the  comparison  of  performance  in  terms  of  mean-square  error 
and  convergence  time.  These  steps  are  to  be  conducted  for  both  the  common  transversal 
structure  and  the  more  complex  lattice  structure. 
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2)  Investigation  of  recursive  adaptation  schemes.  Both  the  transversal  and  lattice 
structures  referred  to  above  are  all-zero  filters,  and  consequently  are  unconditionally  stable 
(although  the  adaptation  process  can  be  unstable).  In  many  cases,  the  tradeoff  for  this 
unconditional  stability  is  increased  complexity  for  a  given  performance  level  in  comparison 
to  the  performance  offered  by  recursive  filters  which  have  both  poles  and  zeros. 
Consequently,  this  task  has  the  goal  of  developing  algorithms  for  adaptive  recursive  filters. 

3)  Unification  of  the  theory.  Where  possible,  the  analyses  conducted  in  this  project 
are  applied  to  both  DT  and  CT  implementations,  with  numerous  comparisons  of  the  two 
approaches  drawn  along  the  way.  The  result  of  this  unification  is  a  body  of  information 
which  can  be  exploited  by  the  engineer  to  address  the  problem  of  adaptive  filtering  without 
being  forced  into  a  digital  or  an  analog  implementation  at  the  outset. 

Activity 

Funding  from  the  Engineering  Foundation  for  this  project  provides  partial  support 
for  the  principal  investigator  (PI),  Dr.  Frank  Cornett,  at  an  intended  level  of  25%  during 
the  academic  year  and  100%  for  two  months  of  the  summer  of  1991.  The  actual  level  of 
effort  for  the  PI  was  closer  to  40%  during  the  academic  year,  made  possible  by  the  release 
of  the  PI  from  the  teaching  of  one  course  during  both  the  fall  and  spring  semesters.  This 
release  has  been  accomplished  with  matching  funds  provided  by  the  Center  for  Manufactur¬ 
ing  Research  and  Technology  Utilization  at  Tennessee  Technological  University.  The  PI 
utilized  this  time  for  both  conducting  his  own  research  and  for  advising  two  students  who 
conducted  research  leading  to  the  masters  degree.  The  efforts  of  all  three  of  these  individu¬ 
als  are  directly  applicable  to  the  goals  of  this  project. 
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One  student  is  Mr.  David  Reynolds,  who,  having  completed  his  course  work,  devoted 
100%  of  his  time  since  the  inception  of  this  grant  to  research  for  the  completion  of  his 
degree  requirements.  Mr.  Reynolds  completed  the  writing  of  his  thesis,  fulfilled  degree 
requirements  in  March,  and  graduated  officially  in  May  of  1991.  Although  his  research 
applies  to  this  project,  Mr.  Reynolds  was  not  supported  under  this  grant 

The  second  student  under  the  advisement  of  the  PI  is  Mr.  Arunkumar  Rajasekaran, 
who  began  his  masters  program  in  the  fall  of  1990.  Mr.  Rajasekaran  has  completed  most 
of  his  course  work,  and  is  conducting  research  in  the  area  of  recursive  adaptive  filters. 
During  the  summer  of  1991,  he  worked  full  time  at  his  research.  He  is  supported  by  funds 
from  the  Manufacturing  Center  under  the  terms  of  this  grant. 

Progress 

In  order  to  facilitate  the  discussions  to  follow,  a  block  diagram  illustrating  the  so- 
called  direct  form  structure  of  the  continuous-time  adaptive  transversal  filter  is  provided  in 
Figure  1.  The  filter  consists  of  a  sequence  of  analog  blocks,  each  with  transfer  function 
H(s),  feedforward  coefficients  a,,  i=0,l,...,M,  and  feedback  coefficients  bj,  j  =  1,2,...,N.  For 
purposes  of  illustration,  M  =  N  =  3  in  Figure  1,  although  these  parameters  are  unrestricted 
in  practice.  The  process  which  adapts  the  a  and  b  coefficients  in  order  to  meet  the  filter 
performance  criterion  is  not  shown  explicitly.  Those  familiar  with  DT  adaptive  filters  will 
note  that  the  structure  in  Figure  1  is  identical  to  that  for  the  DT  implementation,  except 
that  in  the  DT  structure,  H(s)  is  replaced  by  H(z)  =  z'1,  the  unit  delay  operator.  This 
similarity  in  structure  leads  to  a  number  of  commonalities  in  performance  between  the  CT 
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FIGURE  1 

DIRECT-FORM  STRUCTURE 


and  DT  adaptive  filters.  An  equivalent  similarity  and  set  of  commonalities  exists  between 
the  CT  and  DT  lattice  structures[Ref.  1]. 

Progress  to  date  can  be  stated  in  terms  of  the  three  objectives  delineated  above: 

1)  Comparison  of  discrete-time  and  continuous-time  adaptive  filters.  All  of  the  work 
conducted  within  this  task  has  been  confined  to  transversal  filters,  that  is,  filters  for  which 
the  feedback  coefficients  b4  in  Figure  1  are  all  zero,  and  to  lattice  filters  with  only 
feedforward  signal  propagation.  Investigation  of  recursive  filters,  that  is,  filters  with 
feedback,  was  carried  out  as  a  separate  task,  described  below. 

Analyses  were  conducted  in  the  context  of  interference  reduction  or  signal 
equalization,  both  of  which  are  common  signal  processing  operations  in  the  field  of 
telecommunications.  In  these  contexts,  the  input  signal  x(t)  shown  in  Figure  1  consists  of 
a  desired  signal  (possibly  distorted  by  an  unknown  filter)  plus  noise  and/or  interference. 
The  role  of  the  adaptive  filter  is  to  minimize  the  power  (mean-square)  in  the  error  signal 
e(t)  =  y(t)  -  d(t),  where  d(t)  is  the  desired  signal,  undistorted  and  noise  free,  and  y(t)  is  the 
filter  output 

For  the  direct  form  adaptive  structure,  the  following  progress  has  been  made:  an 
adaptation  algorithm  has  been  derived  which  continuously  updates  the  feedforward 
coefficients  aj  so  as  to  minimize  the  mean-square  error  (MSE).  Two  critical  performance 
characteristics  in  the  adaptation  process  are  the  time  required  for  the  filter  coefficients  to 
converge  to  their  proper  values,  and  the  value  of  the  minimum  MSE  (MMSE)  which  can 
be  achieved  after  convergence  is  complete.  Analysis  of  both  characteristics  has  been 
completed,  showing  that  the  performance  of  the  CT  adaptive  filter  is  essentially  identical 
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to  that  of  the  DT  filter.  This  result  is  critical,  indicating  that  the  design  engineer  need  not 
give  up  performance  in  order  to  obtain  the  bandwidth  and  complexity  advantages  of  the 
analog  implementation. 

Sensitivity  of  the  MMSE  to  errors  in  the  bandwidth  of  the  basic  analog  block  H(s) 
has  been  determined.  The  effects  of  such  errors,  which  are  inherent  in  an  analog  system, 
are  substantially  mitigated  by  the  adaptation  process  itself  in  an  adaptive  structure.  The 
adaptation  simply  causes  the  filter  coefficients  to  be  modified  in  such  a  way  as  to 
accommodate  the  bandwidth  errors,  in  addition  to  minimizing  the  MSE.  Consequently,  if 
the  parameters  of  H(s)  are  controlled  to  reasonable  tolerances,  as  would  be  expected, 
inaccuracies  do  not  result  in  significant  system  performance  degradation. 

For  the  CT  lattice  structure,  both  a  normalized  and  an  unnormalized  adaptation 
algorithm  have  been  derived.  Convergence  and  MMSE  characteristics  have  been 
determined,  once  again  indicating  that  the  CT  structure  yields  performance  identical  to  the 
DT  implementation. 

For  both  the  direct  form  and  lattice  structures,  analyses  have  been  completed  to 
determine  the  effect  of  the  choice  for  the  basic  analog  block  H(s)  used  in  the  adaptive  filter. 
Both  one-pole  lowpass  and  one-pole  allpass  blocks  have  been  considered,  and  the  MMSE 
has  been  shown  to  be  identical  for  both  block  types.  However,  adaptive  structures  made 
up  of  lowpass  blocks  exhibit  in  many  situations  a  very  large  spread  between  the  maximum 
and  minimum  eigenvalues  of  the  signal  autocorrelation  matrix,  which  results  in  very  slow 
convergence.  By  comparison,  the  eigenvalue  spread  for  structures  incorporating  allpass 
blocks  may  be  smaller  by  several  orders  of  magnitude,  indicating  that  the  allpass  block  is 


103 


preferable  in  order  to  obtain  good  convergence  characteristics.  A  second  reason  for 
preferring  the  allpass  block  over  the  lowpass  block  is  that  the  coefficients  for  the  filter 
made  up  of  allpass  blocks  have  a  much  smaller  minimum-to-maximum  range.  Since  a  large 
dynamic  range  is  difficult  to  achieve  in  analog  multipliers,  this  result  is  very  important  in 
terms  of  a  practical  implementation. 

In  addition  to  allpass  and  lowpass  characteristics,  the  analog  block  H(s)  can  be 
chosen  to  be  an  integrator,  for  which  H(s)  =  1/s,  or  a  differentiator,  for  which  H(s)  =  s. 
These  two  choices  lead  to  a  simplification  in  the  analysis  of  the  adaptive  structure,  and  also 
result  in  no  interaction  of  the  adaptive  filter  coefficients,  as  occurs  with  either  an  allpass  or 
lowpass  basic  analog  block.  There  are  numerous  practical  difficulties  associated  with 
integrators  and  differentiators,  in  particular,  differentiators  exhibit  pioblems  with  high- 
frequency  noise,  and  integrators  are  plagued  by  nonzero  offset  voltages.  Despite  these 
difficulties,  it  was  felt  that  the  potential  for  simplification  warranted  an  investigation  of  their 
use.  In  the  investigation  conducted,  which  was  for  a  particular  application  of  equalization, 
much  the  same  set  of  problems  was  encountered  as  for  the  lowpass  analog  block,  leading 
to  the  conclusion  that,  among  the  four  types  of  basic  analog  block  types  considered,  the 
allpass  structure  provides  the  best  performance.  It  is  also  practical  to  implement,  consisting 
of  a  simple  RC  network  and  an  active  amplifier  with  a  gain  of  2. 

This  preference  for  the  allpass  block  carries  over  to  the  lattice  structure  as  well.  An 
important  benefit  of  the  lattice  filter  is  the  orthogonalization  of  the  backward  error  residual 
signals.  This  feature  occurs  naturally  in  the  DT  lattice  filter,  but  is  obtained  in  the  CT  filter 
only  if  the  analog  blocks  are  allpass. 
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Numerous  computer  simulations  have  been  completed  for  both  analog  and  digital 
implementations  of  both  direct  form  and  lattice  structures.  These  simulations  corroborate 
the  theoretical  analysis  very  well,  and  bear  out  the  performance  equivalence  between  digital 
and  analog  adaptive  filters. 

During  the  summer  of  1991,  the  results  achieved  within  this  task  were  incorporated 
into  two  papers  which  are  planned  for  submission  for  puolication  in  the  Transactions  on 
Circuits  and  Systems  of  the  IEEE.  The  first  paper  delineates  results  related  to  the  direct 
form  structure,  whereas  the  second  provides  results  for  the  lattice  structure. 

2)  Investigation  of  recursive  adaptation  schemes.  It  is  well  known  that  in  most 
filtering  applications,  recursive  filters,  with  both  poles  and  zeros,  provide  better  performance 
than  nonrecursive,  or  all-zero  filters  for  a  given  filter  complexity,  where  complexity  can  be 
measured  as  the  total  number  of  feedforward  and  feedback  coefficients.  However,  a 
number  of  problems  arise  with  recursive  adaptive  filters.  First  is  the  problem  of  filter 
stability.  Although  the  adaptation  process  for  nonrecursive  filters  can  be  unstable,  instability 
is  seldom  a  problem  since  it  can  be  avoided  with  small  adaptation  gain.  Instability  of  the 
filter  itself  is  impossible.  With  recursive  filters,  there  is  a  potential  for  instability  in  both 
the  adaptation  process  and  the  fiUer.  Secondly,  the  error  surface  for  the  recursive  filter 
exhibits  local  mimima,  so  that  convergence  of  the  filter  coefficients  to  attain  the  global 
minimum  of  the  MSE  is  !u  general  not  assured.  Thirdly,  the  adaptation  process  for 
recursive  filters  is  more  complicated,  for  a  given  complexity,  than  for  nonrecursive  filters. 

In  the  late  1970’s,  several  papers[2,3,4,5]  appeared  describinp  the  so-called 
hyperstable  adaptive  recursive  filter  (HARF).  This  technique  was  developed  in  the  context 
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of  identifying  the  pole-zero  coefficients  of  an  unknown  discrete-time  plant  However,  by 
recasting  the  problem,  HARF  can  be  applied  in  equalization  or  interference  reduction 
applications.  The  HARF  algorithm  solves  the  first  two  problems  described  above,  making 
it  an  attractive  alternative  to  the  nonrecursive  adaptive  scheme.  A  significant  reduction  in 
complexity  was  proposed[4,5]  as  the  simplified  HARF,  or  SHARP.  The  cost  of  thisTeduced 
complexity  is  that  the  SHARF  cannot  be  shown  to  be  stable. 

In  addition  to  the  above  mentioned  shortcomings  of  HARF  and  SHARF,  both 
algorithms  share  a  second  serious  drawback.  If  the  unknown  plant  has  transfer  function 
B(z)/A(z),  then  a  smoothing  filter  C(z)  must  be  employed  in  the  HARF  adaptation  process, 
where  the  transfer  function  C(z)/A(z)  has  phase  between  -90°  and  90°  for  all  values  of  z  on 
the  unit  circle.  Thus,  the  proper  choice  for  C(z)  requires  some  knowledge  of  A(z),  which 
is  often  unavailable. 

As  part  of  this  project,  a  CT  HARF  algorithm  has  been  developed.  This  algorithm 
shows  great  promise,  since  it  naturally  has  the  reduced  complexity  of  the  DT  SHARF 
technique.  More  importantly,  stability  is  not  sacrificed  with  the  CT  HARF  filter;  a  proof 
demonstrating  convergence  has  been  developed.  A  limited  number  of  simulations  have  been 
conducted. 

Unfortunately,  the  CT  HARF  shares  with  the  DT  algorithms  the  drawback  of 
requiring  a  smoothing  filter  C(s)  for  which  the  phase  of  C(s)/A(s)  must  lie  between  -90°  and 
90°  for  all  s  on  the  jw  axis.  Consequently,  modified  HARF  algorithms  have  been  developed, 
for  both  CT  and  DT  applications,  which  obviate  the  need  for  knowledge  of  A(z)  or  A(s). 
Both  adaptation  algorithms  have  been  shown  to  be  stable,  meaning  that  the  error  between 
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the  desired  signal  and  the  adaptive  filter  output  approaches  zero  as  time  progresses. 
However,  because  the  filter  output  is  not  computed  in  a  truly  recursive  manner  in  the 
modified  algorithms,  it  is  not  for  certain  that  the  filter  which  results  from  the  adaptation  is 
stable.  Thus,  further  analysis  and  simulation  will  be  conducted  to  evaluate  these  algorithms. 

Another  approach,  derived  from  the  CHARF(6]  algorithm,  has  been  investigated  to 
a  small  extent  It  appears  that  with  this  approach  the  smoothing  filter  C(s)  can  either  be 
obtained  adaptively  or  eliminated  altogether.  Further  analyses  and  simulations  are  planned 
to  determine  the  effectiveness  of  this  algorithm. 

If  these  efforts  toward  an  effective  continuous-time  recursive  adaptive  filter  prove 
successful,  it  is  felt  that  the  resultant  algorithms  will  constitute  a  significant  contribution  to 
the  field  of  adaptive  filtering. 

3)  Unification  of  the  theory.  As  can  be  seen  from  the  above  discussions,  every  effort 
has  been  made  in  this  project  to  carry  forward  the  treatment  of  both  analog  and  digital 
implementations  of  adaptive  filters  in  parallel.  The  result  of  such  a  treatment  is  a  body  of 
knowledge  relating  to  the  actual  problem  of  adaptive  filtering,  without  predisposition  toward 
the  form  of  implementation. 
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Abstract 

Fault  simulation  constitutes  an  indispensable  tool  in  ensuring  the  correctness  and  quality  of  manufactured  digital  designs. 
TVaditional  uniprocessor  based  algorithms  for  fault  simulation  have  been  observed  to  be  agonizingly  slow  for  today’s  large 
and  complex  digital  designs.  More  recently,  the  literature  presents  an  approach,  introduced  by  a  few  researchers,  wherein 
the  fault  set  was  partitioned  and  the  digital  design  fault  simulated  for  each  of  the  fault  subset  on  separate  processors  of 
a  parallel  processor  system.  This  approach  is  limited  in  that  it  continues  to  utilize  the  traditional  uniprocessor-based 
algorithm  and  the  results  indicate  that  the  performance  results  are  not  encouraging.  This  paper  introduces,  perhaps 
for  the  first  time,  a  novel,  distributed,  circuit-partitioning  based  algorithm  for  fault  simulation  of  both  combinational 
and  sequential  digital  designs.  In  this  approach,  referred  to  as  NODIFS  in  this  paper,  every  component  in  the  circuit 
is  modeled  as  an  asynchronous,  concurrent,  entity  that  is  fault  simulated  as  soon  as  appropriate  signal  transitions  and 
fault  lists  are  asserted  at  its  input  ports.  The  circuit  is  partitioned  such  that  components  of  every  partition  are  allocated 
to  a  unique  processor  of  the  parallel  processor  system.  Consequently,  a  number  of  components  may  be  concurrently 
fault  simulated  on  multiple  processors  in  NODIFS  implying  significant  increase  in  throughput.  The  advantages  of  this 
approach  are  -  (i)  potentially  very  high  throughput  in  contrast  to  the  traditional  appraoch,  (ii)  the  algorithm  is  novel 
in  that  the  overall  task  of  decision-making  i.e.,  fault  simulation  of  the  circuit,  is  distributed  into  a  number  of  natural, 
independent,  and  concurrent  entities  that  execute  asynchronously  to  utilize  maximal  parallelism,  and  (iii)  the  algorithm 
is  scalable.  The  notion  of  performance  scalability  implies  that  where  the  problem  size  i.e.,  the  size  of  the  circuit  being 
fault  simulated,  is  doubled,  and  the  number  of  processors  of  the  parallel  processor  system  is  simultaneously  increased 
by  a  factor  of  2,  the  increase  in  the  required  CPU  time  is  not  appreciable.  Equivalently,  the  level  of  performance  is 
held  at  a  uniform  value  of  less  than  2.0.  For  a  given  circuit  size,  S,  number  of  processors,  P,  and  speedup  factor  /j, 
the  performance  level  is  defined  in  this  paper  as  the  ratio  of  the  speedup  factor  f?  to  /,,  where  /2  corresponds  to 
the  circuit  size  and  number  of  processors  doubled  to  2xS  and  2xP  respectively.  The  speedup  factor  is  defined  as  the 
ratio  of  the  CPU  times  required  for  fault  simulation  in  the  uniprocessor  to  that  for  the  distributed  approach.  Given 
that  the  ratio  of  the  problem  size  to  the  number  of  computational  engines  may  be  effectively  maintained  at  a  uniform 
value  for  any  arbitrary  yet  pragmatic  problem  size,  the  possibility  of  continued  increase  in  the  overall  performance  of 
the  fault  simulation  algorithm,  as  opposed  to  a  uniprocessor,  is  implied.  Additionally,  the  NODIFS  algorithm  is  a 
natural  choice  for  fault  simulation  of  behavior-level  circuits  -  an  eventual  reality,  wherein  the  ratio  of  the  computational 
to  communication  load  for  the  behavior  models  may  approach  a  significantly  large  value.  This  paper  also  reports  on 
an  implementation  of  NODIFS  on  both  the  ARMSTRONG  parallel  processor  system  at  Brown  University  and  the 
performance  results  indicate  significant  increase  in  the  speedup  for  a  few  representative  example  digital  designs. 
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Dr.  Charles  Fr{Sman 
Engineering  Foundation 
Suite  303 

345  East  47tfc  Street 
New  York.  NY  10017 


Dear  Dr.  FrWman, 

Please  find  enclosed  a  report  that  details  the  results  of  the  second  and  last  phase  of  the  research  conducted 
under  the  grant  EF-RI-B-90-22.  Thus,  I  have  completed  the  goals  as  stated  in  my  letter  to  you  dated 
October  1, 1992.  I  am  very  grateful  to  the  Engineering  Foundation  for  permitting  me  the  no-cost  extension 
to  the  original  research  proposal  for  1991-1992.  I  believe  that  this  second  phase  of  the  research  will  be 
extremely  valuable  to  the  industry  since  most  businesses'lack  specialized  parallel  processors  but  possess 
large  numbers  of  workstations.  I  am  already  in  the  process  of  using  the  contents  of  the  enclosed  draft  to 
prepare  a  journal  quality  paper  for  possible  publication  in  the  IEEE  Computer.  I  hope  to  finish  writing  the 
paper  by  August  1992  and  will  be  delighted  to  acknowledge  the  support  of  this  grant.  I  will  also  mail  you  a 
copy  of  the  final  paper  at  that  time.  I  will  continue  to  pursue  the  spirit  of  this  research,  i.e.  fault  simulate 
even  larger  circuits.  The  first  paper  from  this  grant  is  still  under  review  by  the  International  Journal  of 
Electronic  Testing  (JETTA)  and  a  conference  version  of  the  paper  has  been  published  in  the  proceedings 
of  the  IEEE  Systems,  Man,  and  Cybernetics  Conference,  October  1991. 

Thank  you  very  much, 

Sincerely, 
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as  a  Loosely-Coupled  Parallel  Processor  System 


Kenneth  S.  Lyons  Sumit  Ghosh 
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Rhode  Island  02912 

The  ARMSTRONG  system  is  known  as  a  loosely-coupled  parallel  processor,  each  node  having  its  own 
memory  (local  memory)  and  a  means  of  communication  with  other  nodes.  In  the  ARMSTRONG  case,  each 
node  is  physically  connected  to  up  to  eight  other  nodes,  but  can  communicate  to  any  node  through  message 
routing.  The  physical  connections  between  the  ARMSTRONG  nodes  are  designed  as  patchcords  that  can 
be  configured  by  hand  to  the  desired  network  topology.  This  is  a  great  asset  in  optimizing  the  network.  It 
is  easy  to  see  the  similarities  between  the  ARMSTRONG  parallel  processor  and  a  network  of  computers- 
in  our  case  a  network  of  Sun  workstations.  Each  workstation  has  its  own  memory  and  can  communicate 
with  any  other  workstation.  One  extremely  big  advantage  is  the  fact  that  each  Sun  workstation  is  far  more 
powerful  than  an  ARMSTRONG  node.  It  is  difficult  to  say  exactly  how  much  more  powerful,  but  it  should 
be  noted  that  the  Suns  are  running  a  68030  or  68040  at  xxMHz  with  true  coprocessor  interface,  while  an  A 

In  practice  only  small  simulations  could  be  run  on  the  ARMSTRONG,  as  the  nodes  quickly  ran  out  of 
memory.  At  this  point  it  became  necessary  to  lay  down  the  foundations  for  using  Sun  workstations  as  a 
loosely-coupled  network  of  parallel  processors.  The  project  began  with  modifying  existing  ARMSTRONG 
simulations  so  they  could  be  run  on  the  Suns.  This  was  fairly  easy  as  the  inter-process  comunications  on 
the  ARMSTRONG  system  strongly  resemble  internet  sockets. 

The  code  was  written  as  separate  functions  that  emulate  ARMSTRONG  built-in  functions.  We  began  with 
the  ARMSTRONG  ns_bind  call,  syntax:  ns_bind(processjiame).  This  is  called  by  all  processes  to  assign  a 
symbolic  name  to  themselves.  The  nodes  cannot  receive  specific  communications  without  first  assigning 
themselves  a  name,  but  they  can  receive  datagrams.  To  emulate  this  function  on  the  Suns,  it  is  necessary 
to  save  the  pertinent  information  about  the  calling  process  in  a  file.  The  file  is  called  "boundnames”  and 
consists  of  a  line  for  each  process  that  has  called  ns_bind.  On  this  line  is  the  name  of  the  calling  process 
(the  argument  of  the  ns_bind  call),  followed  by  the  name  of  the  Sun  workstation  it  is  running  on,  followed  by 
the  internet  address  of  the  workstation,  followed  by  a  unique  32  bit  number  that  identifies  the  workstation’s 
CPU,  followed  by  the  local  process  number  of  the  calling  process.  Not  all  of  this  information  is  necessary 
for  full  compatibility,  but  is  provided  for  future  flexibility 


/*  Function  ns_bind(process_name) 

/*  assigns  symbolic  name  process_name  to  calling  process 
/*  called  by:  anyone 
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*/ 

*/ 

*/ 


*/ 

*/ 

*/ 


/*  calls  :  establish_connection()  to  create  a  socket  to  listen  on 
/*  (various  system  calls) 

/*  files  :  boundnames 

int  ns_bind(char  name[MAXNAME] ) 

{ 

struct  hostent  *hp; 
char  hostname [MAXNAME] ; 
int  hn,nlen=sizeof (hostname) ; 
char  hostaddress [20] ; 
long  host id; 
int  process; 

FILE  *fp; 

if  ( (hn=gethostname (hostname, nlen) )  <0)  {  /*  get  name,  of  workstation  */ 

perror( "Couldn't  find  host  name!"); 
return (hn) ; 

> 

hp=getho st byname (hostname) ;  /*  get  internet  address  */ 
if  (hp==NULL)  { 
retum(-l) ; 

> 

hostid=gethostid() ;  /*  get  unique  32  bit  CPU  number  */ 

process=getpid() ;  /*  get  local  process  number  */ 

bcopy  (hp->h_addr, hostaddress, hp->h_length) ; 

while  ((fp=fopen("boundnames","a"))  ==NULL) ; 

/♦keep  trying  file  in  case  someone  else  is  accessing  it  at  the  time  ♦/ 

fprintf  (fp,"Xs  Xs  '/,s  V,u  y,u\n" .name, hostname, inet_ntoa(hostaddress)  ,hostid, process)  ; 
fclose(fp) ; 

establish. connectionO ;  /*  create  a  socket  to  listen  on  */ 
return (0) ; 

> 


Note  that  this  function  makes  a  call  to  establish.connection,  which  has  no  counterpart  on  the  ARMSTRONG 
system.  This  function  is  necessary  to  create  a  socket  on  the  Suns,  and  is  designed  to  return  the  socket 
identifier  in  the  global  variable  globalsocket.  It  is  necessary  to  use  a  global  variable  for  this  purpose, 
as  the  ARMSTRONG  ipc.accept  call  does  not  take  an  argument  for  the  socket  to  accept  on  (since  the 
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ARMSTRONG  does  not  need  sockets).  The  code  for  establish.connection  is  as  follows: 


/*  Function  establish_connection()  */ 
/*  creates  a  socket  to  listen  for  connect  requests  */ 
/*  called  by:  ns_bind  */ 
/*  calls  :  (various  system  calls)  */ 


int  establish_connection(void) 

{ 

struct  sockaddr_in  name; 

/*  create  a  stream  socket  (reliable  communication)  of  type  internet  */ 
if  ( (globalsocketasocket (AF_INET , SQCK_STREAM, 0) ) <0)  { 
perror(" - >  Error  creating  socket") ; 

exit(l);  /*  could  not  open  socket,  no  point  to  continue  execution  */ 

> 

/*  initialize  address  structure  */ 
name . sin_f amily=AF_INET ; 

name.sin_addr.s_addr=INADDR_ANY;  /*  wildcard  address  */ 
name . sin_port=htons (PORT) ;  /*  PORT  should  be  defined  earlier  as 

/*  unused  port 

/*  bind  address  to  socket  */ 
if  (bind (globalsocket , (struct  sockaddr  *)&name,sizeof  name)  <0)  { 
perror(" - >  Error  binding  socket"); 

exit (2);  /*  could  not  bind  address,  no  point  to  continue  execution  */ 

> 

listen (globalsocket ,5) ;  /*  queue  up  to  five  connect  requests  */ 
return(O) ;  /*  exit  ok  */ 

> 


an  */ 
*/ 


Now  that  we  have  shown  how  to  bind  a  name  to  a  process  and  create  a  socket  to  receive  connections, 
we  need  to  add  functions  to  perform  the  accept,  connect,  send  and  receive  routines.  The  ipc.connect  call 
is  fairly  straightforward  on  the  ARMSTRONG,  taking  the  name  of  the  process  to  connect  to  as  its  first 
argument  and  returning  the  connection  identifier  in  its  second  argument.  In  the  socket  analogy,  the  connect 
call  requires  the  internet  address  to  connect  to  and  returns  the  connection  identifier.  The  connect  request 
must  also  be  given  a  socket  descriptor  to  perform  the  connection  on.  This  is  done  by  creating  a  new 
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socket  each  time  the  ipc.connect  function  is  called.  The  simplest  form  of  this  function  would  access  the 
"boundnames”  file  to  find  the  name  of  the  Sun  workstation  that  corresponds  to  the  process  name  given  in 
the  ipc.connect  call,  then  make  a  connection  to  that  workstation.  The  code  to  do  this  is: 


*/ 
*/ 
*/ 

*/ 

int  ipc.connect (char  name [MAXNAME] ,int  *cid) 


/*  Function  ipc.connect (process.name, connect ion.id) 

/*  connects  to  process  with,  process.name  and  returns  connection  identifier 
/*  called  by:  anyone 

/*  calls  :  (various  system  calls)  */ 

/*  files  :  boundnames 


char  processname [MAXNAME] ; 
char  hostname [MAXNAME] ; 
char  hostaddress [20] ; 
char  hostid[20] ; 
char  processIlO] ; 
char  c; 

struct  sockaddr.in  sockname; 
struct  hostent  *hp; 
int  i; 

FILE  *fp; 

while  ((fp-f open ("boundnames" , "r") )  ==NULL)  ;/*  keep  trying  to  open  file  */ 

/*  in  case  someone  else  is  using  it  */ 

while  (1)  { 

fscanf  (fp,"%s  7,s  %s  7. s  Xs"  .processname, hostname, hostaddress, hostid, process)  ;  /*  read 
if  ( Istrcmp (processname ,name))  {  /*  do  we  have  a  match?  */ 
fclose(fp) ;  /*  Yes,  close  file  */ 

sockname. sin_family=AF_INET;  /*  initialize  address  structure  */ 
sockname. sin_port=htons (PORT) ; 

hp  =  gethost byname (hostname) ;  /*  get  internet  address  of  host  */ 
bcopy  (hp->h_addr,  & (sockname. sin.addr.s.addr) ,  hp->h_length) ; 

/*  create  a  socket  to  make  this  connection  on  */ 

if  ( (sock=socket (AF.INET , SOCK. STREAM, 0) ) <0)  { 

perror(" - >  Error  creating  socket") ; 

exit(l);  /*  no  point  in  continuing  */ 

> 
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/*  make  the  connection  */ 

if  (connect(sock,&sockname,sizeof (sockname))  <  0)  { 

close(sock) ; 

return(-l);  /*  no  good,  had  an  error  */ 

> 

*cid=sock; 

return (0) ;  /*  everything’s  fine  */ 

} 

/*  check  for  end  of  file  */ 
c=getc(fp);  /*  should  be  newline  char  */ 
c=getc(fp) ; 
if  (c==E0F)  { 
fclose(fp) ; 

return(-l);  /*  process  not  found,  return  with  error  */ 

> 

ungetc(c,fp) ;  /*  not  end  of  file  so  put  back  good  char  */ 

> 

> 


The  ipc-send  and  ipc_recv  calls  are  a  little  more  difficult,  as  they  can  handle  a  variable  length  argument 
list,  terminated  by  the  flag  EOA.  The  best  way  to  handle  these  routines  is  to  write  several  of  them,  each 
one  specifically  for  sending  or  receiving  a  certain  number  of  arguments  which  are  of  a  particular  type.  Of 
course  this  could  grow  to  be  quite  a  nuisance,  but  to  date  most  of  the  ARMSTRONG  simulations  involve 
only  one  or  two  different  send  and  receive  calls.  In  the  future  when  programs  are  written  directly  on  the 
Suns  this  will  not  be  a  problem;  it  is  only  a  problem  when  converting  ARMSTRONG  simulations  to  run 
on  the  Suns.  The  following  example  is  very  specific.  It  shows  how  the  code  would  be  written  to  send  and 
receive  two  integers  at  a  time.  The  ARMSTRONG  calls  would  be: 

ipcjsend(cid,INT(x),INT(y),EOA) 

Where  cid  is  the  connection  identifier  to  send  on,  INT  is  a  macro  that  specifies  x  and  y  as  integers,  and 
EOA  is  the  'end-of-list'  flag.  The  receive  call  would  be: 

ipc_recv(cid,RINT(x),RINT(y),EOA) 

With  similar  parameters  except  the  macro  RINT  which  tells  the  computer  that  it  is  receiving  an  integer.  To 
implement  these  routines  on  the  Sun,  we  would  need  the  following  two  functions  designed  to  handle  two 
and  only  two  integers: 


#define  INT(x)  x  /*  'dummy'  macros  to  make  functions  compatible  */ 

#def ine  RINT(x)  *x 
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/*  Function  ipc_send(cid,INT(x) ,INT(y) ,E0A) 
/*  sends  two  integers  on  connection  cid 


int  ipc_send(int  protocol_id,int  xl,int  x2,int  dummy) 
int  a; 

a=send (protocoled, &xl,sizeof  (xl)  ,0)  ; 
if  (a<0) 

return(a);  /*  error  */ 

a=send (protocol. id, &x2,sizeof (x2) ,0) ; 
if  (a<0) 

retum(a) ;  /*  error  */ 

return (0) ;  /*  good  */ 

> 

/*  Function  ipc_recv(cid,RINT(x) ,RINT(y) ,E0A)  */ 

/*  receives  two  integers  from  connection  cid  */ 

int  ipc_recv(int  protocoled, int  *xl,int  *x2,int  dummy) 

{ 

int  a; 

a«recv (protocoled, xl  ,sizeof  (int)  ,0)  ; 
if  (a<0) 

return(a) ;  /*  error  */ 

a=recv(protocol_id,x2,sizeof (int) ,0) ; 
if  (a<0) 

return (a) ;  /*  error  */ 

retum(O);  /*  good  */ 

> 


These  function  have  been  expanded  to  include  network  debugging  and  are  joined  by  the  functions  n- 
s.unlookup,  ipc-select,  and  ipc.accept  in  the  file  Armstrong.to.Sun.c  which  can  be  inserted  into  any  ARM¬ 
STRONG  simulation  to  make  it  run  on  the  Suns  (after  modifying  the  ipc_send  and  ipc.recv  calls  of  course). 
This  code  is  included  in  Appendix  A. 

As  mentioned  earlier,  the  Sun  workstations  communicate  over  ethernet  in  a  ring  network  through  message 
routing.  It  is  possible  that  this  network  could  become  quite  congested  by  running  large  simulations,  especially 
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if  the  network  is  being  shared  with  other  users.  The  ARMSTRONG  computer  allows  you  to  physically 
construct  the  desired  network  as  a  tree,  hypercube,  or  other  structure,  but  the  Suns  are  basically  limited 
to  the  ring  structure.  To  determine  the  feasibility  of  using  the  Suns  for  parallel  simulation,  a  benchmark 
program  was  written  to  send  packets  to  all  other  nodes  running  the  benchmark  program.  The  code  for  the 
program  and  the  results  of  the  simulation  are  found  in  Appendix  B. 


Appendix  A:  Armstrong_to_Sun.c 

/*  This  is  a  file  of  simple  routines  that  can  be  added  to  any  Armstrong 

*  program  to  make  it  run  on  the  Suns,  with  little  or  no  modification. 

*  Note  that  the  functions  are  designed  to  mimic  ipc  calls  on  the  Armstrong 

*  by  using  sockets  on  the  Suns.  All  of  the  code  is  modular,  so  it  should 

*  be  easy  to  upgrade  it  to  STREAMS  or  other  forms  of  communication. 

*  Note  that  the  function  establish_connection  does  not  have  a  counterpart 

*  on  the  Armstrong  system.  It  is  necessary  to  initialize  a  socket  on  which 

*  to  listen  for  connect  requests,  and  that  is  accomplished  by  this  function. 

*  It  needs  to  be  called  only  once  before  any  ipc_accept  calls  are  made,  and 

*  is  called  automatically  inside  of  an  ns_bind  call  or  an  ipc_accept  call 

*  if  it  has  not  been  called  already.  This  allows  you  to  accept  connections 

*  even  if  you  have  not  assigned  yourself  a  name;  even  though  this  would  be  a 

*  paradox  on  the  Armstrong  system,  it  is  possible  on  the  Suns. 

* 

*  Most  Armstrong  calls  return  a  zero  if  successful,  or  an  error  number  if 

*  unsuccessful.  These  routines  have  been  designed  to  do  the  same,  although 

*  this  limits  the  power  that  Sun  socket  calls  usually  posess. 

* 

*  The  ipc_send  and  ipc_recv  routines  are  very  program  specific,  as  on  the 

*  Armstrong  it  is  possible  to  send  any  number  of  any  type  of  arguments 

*  delimitted  by  an  EOA,  ex: 

* 

*  ipc_send(pid,INT(x) ,STRUCT(y) ,CHAR(a) , STRING (b) ,D0UBLE(z) ,E0A) ; 

*  ipc_recv(pid,RINT(x) ,RSTRUCT(y) ,RCHAR(a) ,RSTRING(b) ,RD0UBLE(z) ,E0A) ; 

* 

*  The  corresponding  calls  on  the  Suns  do  not  allow  this  kind  of  flexibility, 

*  so  you  may  have  to  write  your  own  routines,  one  for  each  different  type 

*  of  send  and  receive. 

* 

*  Example  calls : 

* 

*  establish_connection ()  returns  a  zero  if  successful 

*  ns_bind(pname)  assign  symbolic  name,  returns  zero  if  successful 

*  ns_lookup(pname ,&pid)  finds  global  id  for  process  with  symbolic  name  'pnarae' 
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* 


and  returns  information  in  pid  (structure  gpid) ,  returns  zero  if 

*  successful 

*  ns_unlookup(&pid,pname)  finds  name  for  process  with  global  id  contained  in 

*  pid  and  returns  name  in  pname,  returns  zero  if  successful 

*  ipc_accept(&cid)  accepts  a  connection,  puts  descriptor  in  cid,  returns 

*  a  zero  if  successful 

*  ipc_connect (name,&cid)  connects  to  process  with  symbolic  name  'name', puts 

*  descriptor  in  cid,  returns  a  zero  if  successful 

*  ipc.select (num_cids ,cid_list ,&timo,&idx)  checks  num_cids  descriptors 

*  in  cid_list  within  time  specified  by  timo  for  data  to  be  read. 

*  returns  index  into  cid_list  in  idx,  and  returns  a  zero  if 

*  successful 

*  ipc_send(cid,INT(x) ,INT(y) ,E0A)  send  two  integers  on  cid,  returns 

*  a  zero  if  successful 

*  ipc_recv(cid,RINT(x) ,RINT(y) ,EOA)  receive  two  integer  from  cid,  returns 

*  a  zero  if  successful 

* 

* 

* 

♦/ 


Implementation  and  Performance 

The  distributed  fault  simulation  algorithm,  NODIFS,  has  been  implemented  on  a  network  of  22  SUN 
spare  1-f,  spare  2,  and  3/60  workstations,  utilizing  the  techniques  outlined  in  this  paper.  The  nature  of 
performance  of  NODIFS,  i.e.  the  nature  of  the  graphs  for  the  different  representative  examples,  is  observed 
to  be  similar  to  those  obtained  for  Armstrong.  The  primary  strengths  of  the  new  approach  include  the 
significant  ease  of  debugging  and  the  ability  to  simulate  large;  circuits  while  its  only  limitation  is  the  slower 
communication  speed  arising  from  the  use  of  ethernet  and  TCP/IP.  As  an  extension  to  the  original  grant, 
in  the  past  year,  we  have  successfully  proved  our  original  hypothesis,  namely  the  use  of  the  network  of 
SUN  workstation  as  a  large  loosely-coupled  parallel  processor.  We  could  not,  however,  fault  simulate  larger 
circuits  because  we  ran  out  of  financial  resources  for  student  support. 
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Abstract 

Computer  vision  systems  employ  a  sequence  of  algorithms  that  exhibit  different  computational  char¬ 
acteristics.  These  algorithms  require  different  data  decomposition  and  load  balancing  techniques 
for  efficient  parallel  implementations.  This  paper  presents  several  techniques  to  perform  static  and 
dynamic  data  decomposition  for  common  computer  vision  algorithms.  They  exploit  the  distribu¬ 
tion  of  features  as  a  measure  of  load  for  data  decomposition.  In  these  techniques,  the  distribution 
of  features  computed  during  the  parallel  execution  of  the  current  algorithm  allows  an  informed 
partitioning  for  the  next  algorithm  in  the  pipeline,  keeping  overhead  involved  in  estimating  the 
load  small.  Performance  results  obtained  from  a  shared  memory  multiprocessor  implementation 
using 
given 
up  to 


these  techniques  are  presented.  Furthermore,  a  classification  of  common  vision  algorithms  is 
based  on  their  suitability  for  one  or  more  data  decomposition  techniques.  Improvements  of 
4  times  over  the  performance  of  uniform  block  oriented  partitioning  were  obtained. 
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RADAR  CROSS-SECTION  OF  PARALLEL  CIRCULAR  CYLINDERS: 
A  COMPARATIVE  STUDY  OF  DIFFERENT  TECHNIQUES 


Atef  Z.  Elsherbeni 

Electrical  Engineering  Department 
University  of  Mississippi 
University,  MS  38677 


ABSTRACT 

Numerous  applications  require  accurate  analysis  for  predicting  the  radar  cross- 
section  of  parallel  cylinders.  Among  these  applications  is  the  modelling  of 
reflecting  and  absorbing  surfaces  by  an  array  of  circular  cylinders.  The  problem 
of  the  scattering  by  many  cylinders  has  been  under  investigation  for  decades  using 
different  techniques.  In  this  paper,  the  analysis  of  selected  techniques  will  be 
outlined  along  with  numerical  comparisons  for  the  far  and  near  fields  for  an  array 
of  circular  cylinders.  These  techniques  include  an  asymptotic  high  frequency 
method,  a  boundary  value  approach,  an  iterative  scattering  method,  and  a  hybrid 
method  combining  an  eigenfunction 
expansion  method  with  the  method 
of  moments.  The  parallel  cylinders 
are  arbitrarily  located  as  shown  in 
the  figure  and  are  made  of 
conducting,  dielectric  materials  or 
combinations  of  both.  The  validity 
and  accuracy  of  each  technique  are 
demonstrated  by  illustrative 
examples.  The  limitations  of  the 
associated  computer  codes  when 
using  a  personal  IBM-PC  or 
compatible  computer  are  also 
outlined.  The  effects  of  various 
parameters  (cylinder  to  cylinder 
coupling,  angle  of  incidence  of  the 
electromagnetic  wave,  shadowing 
effects,  spacing  between  the  cylinders  and  the  physical  parameters  of  each 
cylinder)  on  the  radar  cross-section  of  the  array  of  cylinders  are  also  investigated. 


This  work  was  supported  by  the  Air  Force-Engineering  Foundation  Grant  number 
RI-B-90-24. 
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ELECTROMAGNETIC  SCATTERING  FROM  AN  ECCENTRIC 
MULTI-LAYERED  CIRCULAR  CYLINDER 

Ahmed  A.  Kishk,  Rajan  P.  Parrikar,  and  Atef  Z.  Elsherbeni 
Department  of  Electrical  Engineering 
University  of  Mississippi 
University,  Mississippi  38677  U.S.A. 

ABSTRACT 


The  electromagnetic  scattering  from  an  eccentric  multi-layered,  homogeneous,  circular  cylinder 
is  derived  rigorously  by  using  a  boundary  value  mode  matching  approach.  The  layers  are  assumed 
to  be  infinite  in  length  and  their  axes  are  parallel  to  each  other  but  their  centers  may  have  arbitrary 
coordinates.  The  composite  arrangement  is  illuminated  by  a  cylindrical  wave  from  a  line  source 
of  infinite  extent  in  a  direction  parallel  to  the  cylinder  axis.  The  problem  is  two  dimensional  and 
the  solution  to  both  types  of  polarization  is  presented.  The  effects  of  various  geometrical  and 
electrical  parameters  on  the  near  and  far  field  patterns  are  examined. 


This  work  was  partially  supported  by  the  National  Science  Foundation,  Grant  No.  ECS-8906807, 
and  the  Air  Force  Engineering  Foundation,  Grant  RI-B-90-24. 
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CONDUCTING  AND  DIELECTRIC  CYLINDERS 


A.  Z.  Elslierbeni*,  G.  Tian*  and  M.  Hamid*’ 

‘Department  of  Electrical  Engineering  "Department  of  Electrical  Engineering 
University  of  Mississippi  University  of  South  Alabama 

University,  MS  38677  Mobile,  A L  36688 

ABSTRACT 

In  this  paper,  attention  is  focused  on  the  problem  of 
scattering  of  a  Gaussian  beam  from  an  array  of  parallel  circular 
cylinders.  The  cylinders  are  made  of  homogeneous,  isotropic 
dielectric  material,  perfectly  conducting  material,  or  combination. 

An  iterative  scattering  technique  is  developed  to  evaluate  both  the 
near  and  far  field  components.  The  iterative  technique  starts  with 
exact  series  representation  of  the  initial  scattered  fields  from  each 
cylinder  due  to  the  incident  Gaussian  beam.  Then,  the  sum  of  the 
initial  scattered  fields  from  all  cylinders  is  considered  to  be  an 
incident  wave  on  each  individual  cylinder.  Additional  scattered 
field  contribution  is  then  obtained  and  successive  scattering  can 
take  place  until  the  boundary  conditions  on  the  surface  of  every 
cylinder  are  satisfied.  Transverse  electric  (TE)  and  transverse 
magnetic  (TM)  polarized  waves  are  considered  and  the  effects  of 
the  incident  beam  parameters  on  the  far  scattered  field  pattern  are 
investigated. 
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2.  Summary  of  Research 


2.1  Plan  of  Research 

As  described  in  the  original  proposal:  The  research  can  be  divided  into  the  following  six  tasks: 

1 .  Develop  new,  closed  loop,  acquisition  architecture 

2.  Determine  performance  of  new  architecture 

3 .  Investigate  feasibility  of  implementing  the  new  architecture 

4.  Extend  theory  to  the  non-coherent  case 

5 .  Repeat  performance  and  implementation  studies  for  the  non-coherent  case 

6 .  Extend  analysis  to  higher  order  loops 

These  tasks  will  be  executed  sequentially  with  the  amount  of  effort,  risk,  and  potential  benefits 
increasing  with  each  task.  It  is  reasonable  to  expect  that  items  1  through  3  can  be  accomplished  in 
the  first  year  of  research.  Tasks  4  and  5  represent  another  major  step  in  the  development  and  will 
be  initiated,  but  may  not  reach  fruition  in  the  first  year.  The  extension  to  higher  order  loops  will  be 
considered,  but  it  is  difficult  to  be  optimistic  about  developing  meaningful  results  in  a  short  period 
of  time.  The  analysis  of  other  types  of  second,  and  higher,  order  loops  has  proven  to  be  a 
formidable  challenge. 

2.2  Project  Status: 

Tasks  1  and  3  have  been  completed,  as  planned.  The  new  architecture  has  been  developed,  and  the 
feasibility  of  implementation  investigated.  Substantial  progress  has  also  been  made  on  Task  2,  and 
work  is  continuing  in  this  area,  further  details  are  contained  in  Section  4  of  this  report  As 
anticipated  in  the  original  proposal,  Tasks  4  through  6  were  initiated,  but  have  not  yielded  any 
significant  results  to  date.  Work  will  also  continue  on  these  tasks. 
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